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Determination of Strain Distribution 


by the Photo-Grid Process 


GIVEN A. BREWER ano ROBERT B. GLASSCO 
Lockheed Aircraft Corporation* 


ABSTRACT 


The purpose of this paper is to describe the photo-grid process, 
a new method of measuring local deformations. The process 
consists of photographing a pattern of lines upon a sensitized 
specimen which is then tested and the amount of deformation de- 
termined by measuring the new spacing of the lines. Gage 
lengths as small as 0.01 inch are used. By this process the strain 
distribution in tensile specimens of aluminum alloys has been 
plotted; also the strains produced by several forming operations 
have been determined. The maximum local elongation measured 
in tensile test specimens is much higher than the average in two 
inches and is indicative of the elongations obtained in forming 
operations of simple bending or stretching. In deep drawn shells 
even higher strains were measured. The phenomenon of progres- 
sive strain hardening with reference to the strain distribution 
found in the tensile test specimens is discussed. 


INTRODUCTION 


TUDIES OF elongation and ductility of metals and 
S alloys have in the past been confined to results 
measured over comparatively long intervals such as 
the two-inch and eight-inch gage lengths customarily 
used in tensile strength tests. It has been recognized 
that within such gage lengths the elongation may be 
distributed in a non-uniform manner, and that conse- 
quently, the measured overall percentage elongation 
will vary with the gage length used. Such an analysis 
is of prime importance in approaching all types of form- 
ing problems in which the failure of the metal at some 
critical section is the limiting factor. Yet a survey of 
the existing literature reveals that the subject has re- 
ceived little attention. 

In one of the publications' Crane merely alludes to 
the fact that the local elongation at the point of failure 
is much greater than the ‘‘general” elongation in a two- 
inch gage length. In an earlier paper? Krivobok, Lin- 
coln, and Patterson treat this subject more fully; their 
studies were confined to stainless steels of pseudo-aus- 
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Fic. 1. Distribution of elongation, measured every 
quarter inch, over two-inch gage lengths for two austenitic 
stainless steels, both cold-rolled to about 165,000 pounds 
per square inch.? 


tenitic character and are shown in the form of curves of 
distribution of elongation over standard two-inch gage 
lengths (Fig. 1). The effect of strain hardening prop- 
erties on the distribution of local elongations is brought 
out in this reference, and will be given further attention 
later on in the present paper. Elam* is primarily con- 
cerned with the phenomenon whereby the elongation in 
some materials during the tensile test is assumed to oc- 
cur by what may be termed “‘small jumps’’ rather than 
uniformly. True, reference is made to the non-uni- 
formity of strain distribution throughout the specimen 
under test, but because the method of observation was 
quite cumbersome, detailed data were not obtained. 
Templin and Sturm‘ show the variation in overall elon- 
gation when the gage lengths are increased from one 
inch to ten inches. Shorter gage lengths were not 
studied and, therefore, their studies do not contain in- 
formation pertaining to what may be called ‘“‘local dis- 
tribution” of elongation in tensile testing. 

It seems clear that the limitation on previous work 
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has been partly caused by the lack of a suitable tech- 
nique to investigate the local elongations. Experiments 
which may be made to determine the local elongation 
by measuring the reduced cross-sectional area show only 
part of the picture and are limited in application. Even 
the recently developed electric strain gage is not 
adapted to measurements over extremely small inter- 
vals. 

It is, therefore, the purpose of the present paper to 
describe a new method, the photo-grid process, now 
being used at the Lockheed Aircraft Corporation for 
measuring local deformations, and for studying the ulti- 
mate significance thereof in the technical study of metal 


forming. 


TECHNIQUE DEVELOPED IN PRESENT INVESTIGATION 


The experimental technique which is being used in 
the present investigation was developed by Given A. 
Brewer of the Lockheed Aircraft Corporation. Briefly, 
the technique consists of photographing a pattern of 
lines upon a sensitized specimen which is then tested and 
the amount of deformation determined by measuring 
the new spacing of the lines. 

In order to produce the photograph on the specimen, 
a high-contrast negative is first made from the original 
source of the pattern or grid. A very satisfactory 
source for a rectangular coordinate grid is a photoen- 
graver’s glass screen with parallel diamond-ruled lines 
0.01 inch apart. A tensile coupon prepared with this 
pattern is shown in Fig. 2. Other sources such as rec- 
tangular and polar graph paper have been used satis- 
factorily. Exdmples of the use of these two grid 
sources are shown in Fig. 3 and Fig. 4, respectively. 

Following the preparation of the negative, the surface 
of the test specimen is sensitized, and on the sensitized 
surface a contact print of the grid is made. A detailed 
description of this printing process is given in the Ap- 
pendix. 

For examination of the tested specimen, a cathetome- 
ter or measuring microscope is generally used. Occa- 
sionally, it is desirable to make a photo-enlargement 
of the deformed grid and to measure the strain with a 
scale and dividers. This applies to cases in which the 
specimen will not fit in the cathetometer. In these 
latter cases an undeformed ‘‘conttol”’ area is photo- 
graphed in the same picture with the specimen in order 
to establish the scale. 

The special advantages of the technique described 
above are as follows: (1) the printing process eliminates 
all scratches and gage marks which might give stress 
concentrations; (2) the line pattern makes it possible 
to determine the unit deformation of a very small 
length; (3) the use of the engraver’s screen and the 
methods of magnification give good accuracy; (4) two- 
dimensional deformation is shown which is important 
in many cases; (5) the method can be applied to any 
odd-shaped specimen as long as a contact print can be 
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Fic. 2. Photographic enlargement of a tested tensile coupon 
made of 24S-T Alclad sheet, 0.064 inch in thickness, with a 
0.01-inch rectangular photo-grid on one surface. 
made on it before it is deformed; and (6) relatively 
small expense is incurred, whether for large, small, or 
odd-shaped specimens. 
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Fic. 3. Photographic enlargement of a tested tensile coupon 
made of 24S-T Alclad sheet, 0.064 inch in thickness, with a 
0.05-inch rectangular photo-grid on one surface. 
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RESULTS OF TESTS USING THE PHOTO-GRID 


Tensile Tests 


The first application of the photo-grid was on stand- 
ard sheet metal tensile coupons. Fig. 2, already men- 
tioned, shows a typical aluminum alloy tensile speci- 
men on which was imprinted a grid made up from the 
photo-engraver’s screen having lines spaced 0.01 inch 
apart. The same illustration shows an enlargement of 
the section of the specimen adjacent to the rupture. 
It is to be noted that the photographic enlargement 
shows very clearly the details in the area adjacent to 
rupture. Fig. 3 shows a specimen on which was used 
a grid produced with graph paper having 0.05-inch 
squares. By measuring the local unit elongation or 





80 T T if 








4ob/-—-— + 3 











LOCAL UNIT ELONGATION - per cent 


Average over 2° 


oe oe 


POSITION ON CENTER LINE OF SPECIMEN 











/ division = £ inch 


Fic. 5. Distribution of elongation, measured for 0.02- 
inch gage lengths, on the 24S-T Alclad specimen shown 
in Fig. 2. 


strain along the length of these specimens, as already 
described, a strain distribution curve was plotted. In 
Figs. 5 and 6 are shown the curves plotted for the speci- 
mens illustrated in Figs. 2 and 3. Measurements for 
these curves were made every 0.02 inch and 0.05 inch, 
respectively. 

Each point plotted for the above strain distribution 
curves has an abscissa corresponding to the position on 
the undeformed specimen of a small finite increment of 
length as measured from some arbitrary reference 
axis. It has an ordinate equal to the unit elongation 
(strain) of this increment of length as measured from 
the photograph. 

When a curve is drawn through the points, each or- 
dinate then represents unit elongation for a differential 
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Fic. 6. Distribution of elongation, measured for 


0.05-inch gage lengths, on the 24S-T Alclad specimen 
shown in Fig. 3. 
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length of specimen at that location; and the area under 
the curve, being the summation of the elongations of 
each differential length, represents the total elongation 
of the specimen. The area under the curve divided by 
the original length of the specimen then represents the 
average unit elongation, or average strain, over the en- 
tire specimen. This last is the quantity usually deter- 
mined in tensile tests using long gage lengths and is 
commonly expressed as the “‘per cent elongation”’ of the 
material. 

The specimens used as illustrations in this report are 
merely samples taken from the large number of 24S-O 
and 24S-T Alclad aluminum alloys,* specimens which 
have been tested. As might have been expected, there 
has been a large variation in the strain distribution 
curves of the specimens tested; no two are exactly alike. 
However, from the practical standpoint the elongation 
over a short gage length appears important, especially 
in forming small bends, etc. Therefore, data have 
been compiled from many specimens, based on the 
strain distribution curves and showing the elongation 
in '/, inch and 0.01 inch, as well as in the standard 2 


inches. These data show the following ranges of values: 





1/,Inch 0.01 Inch 


Material 2 Inches 
24S-O Alclad 10-22% 25-45% 31-71% 
24S-T Alclad 16-21 25-35% 44-62% 
24S-T Alclad, purchased as 
24S-O and heat treated by 
various methods 10-22% 20-40% 27-65% 





The curves discussed above have shown clearly the 
non-uniform character of the strain distribution 
throughout a specimen when the latter has been pulled 
to failure. What occurs prior to failure has been ob- 
served by taking a series of photographs during the 
pulling of a specimen. These photographs show that 
the pronounced peak in the strain distribution curve 
does not appear until the load is approximately at a 
maximum. The results also show, as might have been 
predicted, that the failure occurs at the location of 
maximum strain, providing that there is no interference 
from extraneous sources. Even though the increments 
of length over which the elongations are measured are 
sufficiently small to show quite accurately the possible 
extent of strain in the alloys under investigation, the 
maximum elongation at the infinitesimal point of 
failure, which might be called the “ultimate strain,’ can- 


not be determined. 


Deep Drawing Tests 


The photo-grid has been successfully applied to tests 
of deep drawing. By using a polar coordinate grid, 
* 24S Alclad has chemical make-up as follows; 4.5% Cu, 
0.6% Mn, 1.5% Mg, remainder Al, and with a pure aluminum 
“O” refers to annealed (soft) temper; ‘‘T’’ refers to 





coating. 


solution heat-treated and aged (hard) temper. 
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measurements were made of both the radial stretching 
and circumferential shrinking occurring over the entire 
surface of the shells. Because of the favorable system 
of combined stresses which occurs in this particular 
forming operation, very high strains are obtainable. 
With 245-O Alclad shells drawn in three stages com- 
pressive strains of 60 per cent and tensile strains of 110 
per cent were obtained without intermediate annealing. 
In Fig. 4 is illustrated one of these specimens—first, the 
unformed blank printed with the polar grid, and second, 
the final shell after three draws. 


Dimpling Tests 


The photo-grid is being used in a study of the dim- 
pling of sheet aluminum parts. This process of flanging 
a rivet hole, when a flush rivet or screw is to be used, 
stretches the rim of the hole. By printing a photo- 
grid of miniature polar coordinates concentrically with 
the hole as illustrated in Fig. 7, the elongations pro- 
duced by both the impact and the squeeze type dim- 
pling tools have been measured. The range of elonga- 
tions obtained in 24S-T Alclad purchased as 24S-O and 
heat treated was found to be from 20 per cent to 32 per 
cent. Maximum permissible values of elongation for 
this forming process are not yet determined, but it is 
interesting to note that these values obtained thus far 
are considerably higher than the average in two inches 
obtained in standard tensile specimens. 


Bending Large Extrusions 


The photo-grid was used for making elongation meas- 
urements on a large’ 24S-T aluminum alloy extrusion 
which had been cracking during a bending operation. 
The measurements showed that this bending operation 
produced a maximum stretching of only 13 per cent in 
the outer fibers. Further measurements showed that 
material of normal ductility (7.e., capable of about 
161/2 per cent elongation over a two-inch gage length 
when pulled as a standard tensile coupon) could suc- 
cessfully withstand a stretching of at least 22 per cent 
when being bent. Consequently, the basic cause of the 
cracking was traced to a portion of the extrusions hav- 
ing subnormal ductility, a result of a work-hardened 
condition along the outer portion which happened to be 
subjected to the greatest stretching. 

In connection with the above study, the measure- 
ments showed how the actual stretching in a bend com- 
pares with the theoretically required elongation calcu- 
lated from the bend radius and the location of the neu- 
tral surface. The calculated elongation for the above 
case was 26 per cent; whereas the actual maximum 
elongation for a 9 degree bend was 13 per cent, and for 
a 30 degree bend, 22 per cent. 

This investigation demonstrated clearly the practical 
importance of a careful study of permissible elongations 
in connection with the design and tooling of curved 
parts. 
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DISCUSSION 


To understand the non-uniformity of strain distribu- 
tion occurring in tensile specimens as illustrated by the 
curves already shown, it must be recognized, that, be- 
cause of the nature of the physical structure of the 
metal, the tensile stress at which yielding begins is not 
exactly the same at every point in the specimen. Fur- 
thermore, the cross-sectional area is not uniform be- 
cause of the difficulty of securing absolutely accurate 
dimensions. Consequently, at some location in the 
specimen the yield resistance, which may be defined as 
the product of area and yield stress, will be a minimum. 
Under the application of a certain load localized yield- 
ing first occurs at this location. Accompanying the 
stretching is a reduction of the cross-sectional area 
which further weakens this section of the specimen. 
Whether the elongation will continue to be concentrated 
at this location, leading ultimately to failure here, or 
whether it will be subsequently distributed in the rest 
of the specimen before failure occurs, depends upon the 
ability of the material to increase its yield resistance at 
the first instant of plastic flow. 

As is well known, practically all metallic alloys, and 
even metals, have the faculty of strain hardening, but of 
widely different magnitude. If the material being 
loaded has an appreciable capacity for strain hardening, 
and if sufficient time is allowed for this phenomenon to 
take place, the reduction of the cross-sectional area at 
the original point of yield (resulting in an increase in 
stress per unit area) is partially compensated by an ac- 
tual increase in the unit resistance of the material. 
The magnitude of the strain hardening may even be 
sufficient to make the section at which the original 
yielding occurred stronger than the rest of the specimen, 
thus temporarily halting yielding of the metal at this 
point. At the next instant the same process will take 





or 





place in some other section where again the yield re- 
sistance is less than the applied load. This cycle is re- 
peated until the rate of strain hardening relative to 
elongation and the rate of area reduction relative to 
elongation are such that the strength of some section no 
longer increases with further stretching. When this 
condition is reached, the local necking down will occur 
and failure is rapidly brought about. Thus, as a result 
of strain hardening, considerable stretching may take 
place at sections which do not break, and a higher 
average elongation is obtained before the specimen fails. 
The shape of the strain distribution curve will obviously 
be different in the two cases. 

The above explanation was offered by Krivobok, 
et al., in their paper already referred to.2. Their studies 
and observations were made on samples of stainless 
steels, in which the degree of physico-chemical insta- 
bility of the alloys, influenced by variations in analysis, 
was found to be the cause of marked differences in local 
strain and its distribution in tensile specimens studied. 
In Fig. 1, for example, it is shown that the low-nickel 
stainless steel ‘A,’ having pronounced work-hardening 
properties, shows a much higher average elongation 
(although the same peak elongation) as the high-nickel 
stainless steel ‘‘B,’’ which does not exhibit the same de- 
gree of strain hardening. 

No phase change is involved in the hardening of 
aluminum alloys, however; and hence physico-chemical 
instability does not enter into the results. It might be 
expected, therefore, that the strain distribution in 
aluminum alloys would differ from that found in cer- 
tain stainless steels. Comparison of Fig. 5 or Fig. 6 
with Fig. 1 indicates that the shape of the strain dis- 
tribution curve for 24S-T Alclad does indeed differ 
from that for stainless steel ““A’’ (the strain-hardening 
steel). It is quite similar, however, to that of stainless 
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two inches shown. 


FIG. 


steel ‘‘B”’ (the non-strain-hardening steel). This agrees 
with the general knowledge that 24S-T exhibits a rela- 
tively slight tendency to strain hardening. 

In a general way, the information presented herein 
leads to the following important conclusions: first, it 
can be seen that for a given material the average unit 
elongation measured after failure is very dependent 
upon the gage length over which it is measured. Second 
it is evident that the customary method of evaluating 
the ductility of a metal, that is, measuring the average 
unit elongation over a standard gage length of say two 
inches, gives a very incomplete picture. It is the shape 
of the strain distribution curve which is of vital impor- 
tance in the study of the forming properties of a mate- 
rial. 

To study further the significance of the shape of the 
strain distribution curve, consider two imaginary met- 
als, ‘‘C’”’ and “D,”’ with the same average strain over 
two inches, but with different maximum local strains. 
Fig. 8A shows the strain distribution curves of these 
metals with the two-inch gage length outlined to show 
the same average over this length for both curves. If 
the length of material being measured is less than two 
inches, say one-quarter of an inch, the average elonga- 
tion of Metal ‘‘C”’ is greater than that of Metal “D,” 
as shown in Fig. 8B. On the other hand, for a gage 
length greater than two inches, say six inches, Metal 
“D” has a higher average elongation as shown in Fig. 
8C. Curves showing average elongation vs. gage 
length for these two metals are given in Fig. 9. 

The fact that the material can actually sustain high 
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Fic.9. Variations of average unit elongation with gage 
length for metals C and D from Fig. 8. 


Distribution of elongation over six-inch gage lengths in two imaginary metals. 
(B) Average unit elongation in one-quarter inch shown. 
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(A) With the average unit elongation in 
(C) Average unit elongation in six inches shown. 


local elongations of magnitude corresponding to the 
peak of the distribution curve, accounts for the fact 
that in sharp bends where elongation is confined to a 
short distance, high strains are found to exist. Such 
strains are often considerably higher than one would be 
led to believe possible on the basis of the allowable elon- 
gations determined over a two-inch gage length. 

It has been found that the average elongation taken 
over a one-quarter-inch gage length on the peak of the 
distribution curve gives a good indication of the relative 
permissible minimum bend radii of various aluminum 
alloy sheets. On the other hand, little or no prediction 
of minimum bend radii can be made on the basis of the 
elongations in two-inch gage lengths. 


PROGRAM OF FUTURE WORK 


It is the desire of the authors that this article be con- 
sidered an invitation to other research engineers to 
adapt the photo-grid technique to the many types of 
experimental work and testing where it can be used to 
advantage. Future work to be conducted using the 
photo-grid should consist of: (1) obtaining strain dis- 
tribution curves on the materials commonly subject to 
forming and other plastic deformation. These curves 
should cover conditions existing throughout the plastic 
range, as well as close to failure; (2) further study of the 
strain hardening phenomenon in relation to the strain 
distribution curves; (3) study of crystal and atomic 
structure as related to the phenomena of strain harden- 
ing and strain distribution; (4) application of the 
knowledge gained above to practical problems of form- 
ing; and (5) application of the photo-grid technique to 
these problems of plastic deformation. 


APPENDIX 
DETAILED DESCRIPTION OF THE PRINTING PROCESS 


In the photo-grid process the printing of the line pat- 
tern onto the specimen is done photographically. The 
main procedures in the printing process are as follows: 
the specimen to be printed is first cleaned with pumice 
and cold water, after which is applied to the surface a 
sensitizing solution consisting of 28 parts by weight of 
water, 4 parts of photoengraving glue, 1 part of am- 
monium bichromate, and !/, part of ammonia water. 
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(This sensitizing solution was a development of Charles 
W. Miller of the Lockheed Aircraft Corporation.) If 
the specimen is a large plate, it is placed on a regular 
photo-engraving plate-whirl consisting of a turn-table 
rotating in warm air, and the solution is poured on. 
If the specimen has an irregular shape, or is too large 
for the plate-whirl, the solution is sprayed on with an 
air gun using 15 pounds per square inch pressure. For 
spraying, 1 ounce of stock solution is dissolved in 4 
ounces of water and 4 ounces of alcohol. 
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PROCESS 


As soon as the surface of the metal has been prepared, 
the negative is placed face to face with the sensitized 
surface and exposed for 5 minutes in a vacuum frame at 
a distance of about 36 inches from a 35 ampere arc 
light. 
with cold water and a solution of Diamond dye (one 


Following the exposure, the specimen is rinsed 


box of Wells and Richardson dye mixed with one pint 
of water) is applied and allowed to dry for 5 minutes, 
after which the specimen is rinsed with cold water and 
again allowed to dry before it is ready for testing. 
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Theoretical Correction for the Lift 
of Elliptic Wings 


ROBERT T. JONESt 


National Advisory Committee for Aeronautics 


ly THE WING section theory the magnitude of the 
circulation, and hence of the lift, is determined by 
the velocity that would be induced near the trailing 
edge of the section in a non-lifting potential flow. In 
three-dimensional flow the problem is complicated by 
the presence of the wake and no simple basic solution 
has been found. Treatment of the problem of a wing of 
finite span has therefore been on the basis of the two- 
dimensional theory, corrected for the effect of the 
wake. * 

The correction commonly applied is that introduced 
by Prandtl and known as the lifting-line theory. The 
downflow induced by the wake is considered, in this 
theory, to reduce the relative normal velocity and hence 
the edge velocity of the wing. It is assumed, however, 
that once the true angle of attack is determined for any 
section of the wing its effect in producing circulation and 
lift is the same as in two-dimensional flow. This as- 
sumption is expressed by the equation 


’ 


Cr = 2r(a — a;) (1) 


where 27 is the slope of the lift curve for the thin wing 
of infinite aspect ratio, a is the angle of attack of the 
section, and a; is the induced angle of downflow. 

Eq. (1) takes into account the effect of the wake in 
diminishing the relative normal velocity of the wing. A 
further correction is indicated by the fact, established in 
hydrodynamic theory, that the surface velocities in- 
duced by a given relative motion of a body in three- 
dimensional flow are generally smaller than those in 
two. For example, the maximum, or “‘edge’’ velocity 
around an infinite cylinder is 1, while that around a 
sphere of the same cross-section is '/2. Similarly, if the 
velocity around the edge of an endless thin plate is 
taken as unity, the corresponding velocity around the 
edge of a circular disc is found to be 2/7, or 0.637. 

In the case of an elliptic disc the velocity at every 
point is reduced by the factor 1/£, where £ is the ratio 
of the semiperimeter of the ellipse to the span. (See 
appendix.) Further investigation shows that the 
chordwise cross-sections of the non-lifting potential 
flow are similar all along the span of the elliptic plate 
and are the same as those obtained from the wing-sec- 
tion theory except for this constant reduction factor. 
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* More recently the problem has been treated on the basis of 
the acceleration potential, which vanishes in the wake.! The 
method is still one of considerable difficulty, however. 


If the velocity distributions of the circulatory flow 
are also assumed to be similar to those given by the 
wing-section theory, the circulation required to satisfy 
the Kutta condition at each section should be reduced 
by the same factor. The corrected formula for the 
lift is then 


Cr = (22/E)(a — a;) (2) 

This correction may be given a physical interpreta- 
tion by considering that a finite wing offers a longer 
edge around which the air may escape (see Fig. 1), 


a 40 > 
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Fic. 1. Correction factor E for the effect of the 
perimeter ratio on the lift: Cz = (297/E)(a@ — a 








and that the air velocities are therefore less in the pro- 
portion that the length of the edge is greater. The rule 
is not exact for plan forms other than the elliptical. 

Since the velocity of the non-lifting potential flow 
is constant all around the edge of the elliptic plate, the 
circulation required will be proportional to the chord at 
each section. The circulation is thus elliptically dis- 
tributed spanwise. Such a distribution, with the chord- 
wise distribution assumed earlier, leads, as in the lifting- 
line theory, to the relation 


a; = C,/rA (3) 


where A is the aspect ratio. 
Substitution of this value into Eq. (2) gives 


Cr = 2raA/(EA + 2) (4) 


Since the chordwise distribution of the circulation in 
three-dimensional flow is assumed similar to that in 
two, and since the similarity is only proved for the non- 
circulatory flow, Eq. (4) must be considered a corree- 
tion of the wing-section theory rather than a solution 
of the three-dimensional problem. The assumption of 
similarity, although its validity is subject to somewhat 
the same limitations in general, appears to be a more 
justifiable one than the assumption of equality made 
in the lifting-line theory. 

The problem of three-dimensional flow around a lift- 
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ing elliptical plate has been treated by Krienes’ using 
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Prandtl’s concept of the acceleration potential. The 
calculations involve the superposition of a series of 
solutions of A*# = 0 in ellipsoidal coordinates (Lamé 
functions), the coefficients of the series being obtained 
from simultaneous equations. Although the overall 
lift depends on only one coefficient, calculations of the 
span loading and the drag would require all the terms 
of the infinite series and was therefore not considered 
practicable. 

The lift coefficients obtained by Krienes’ method are 
shown in Fig. 2. Results obtained by correcting the 
wing-section theory in accordance with Eq. (4) are 
compared with Krienes’, and the agreement is seen to be 
good for ordinary aspect ratios. 

A curve derived from Blenk’s calculation of the vortex 
distribution in a rectangular plate* also agrees closely 
with Eq. (4). 

It is found that the additional correction to the wing- 
section theory accounts for an appreciable fraction of 
the loss in lift that is usually attributed to viscosity. It 
has been difficult to reconcile the magnitude of the in- 
efficiency with the observed dimensions of the wake, 
which in the case of smooth wings, is extremely narrow 
at the trailing edge. The foregoing correction accounts 
for as much as half of this discrepancy in cases of wings 
with sharp trailing edges. (See Fig. 3.) 

According to Munk’s theorem, the induced drag is 
not affected by a displacement of the lifting elements 
in the direction of the chord; hence the relation be- 
tween the lift and the induced drag may be derived on 
the basis of the lifting-line theory. Also, since the pitch- 
ing moment of the elliptic wing involves the same cor- 
rection as the lift, no change of aerodynamic center 
location with aspect ratio is indicated. 

The three-dimensional potential flow around the 
elliptic disc may be used also as the basis for calculating 
the lift of an elliptic wing with varying angle of attack. 
On account of the linearity of the equations, the gen- 
eral problem may be reduced to a determination of the 
growth of lift following a sudden start of the motion with 
the flight velocity V and the normal velocity w = Va.* 

* The infinite force implied by the sudden start leads, in prac- 
tice, to a finite value during the continuous motions to which the 
formulas are applied. 
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Fic. 3. Experimental and theoretical lift-curve slopes. 
Aspect ratio, 6. 


Since the steady or asymptotic value of the lift is 
already known, the most important point is the deter- 
mination of the lift near the start of the motion. 

As a consequence of the Kutta condition, the layer of 
air leaving the wing at the trailing edge must satisfy 
the same boundary condition as the adjacent wing sur- 
face for a short distance downstream. Thus the layer 
of air in contact with the wing retains the motion im- 
parted by the wing and for a short distance after the 
start this layer (the vortex wake) behaves as an im- 
permeable extension of the wing surface.‘ The flow at 
the first instant after the start is thus what might be 
caused by the wing in process of growing wider at the 
rate V while moving with the normal velocity w. The 
starting lift L,. 9 may then be thought of as the reac- 
tion to uniform motion of the wing acting as a body with 
increasing mass: 

Li<o9 = w (dm'/dt) = w (dm'/dc)V (5) 


where m’ 
inertia of the flow and c is the apparent length of the 


is the mass representing the aerodynamic 


chord. In two-dimensional flow 
m' = x (c?/4)p (6) 
so that 
Cig = 7 (7) 


In the case of the elliptic wing the starting process 
is pictured as a change into a similar but slightly wider 
ellipse. For normal aspect ratios this change fulfills 
the required conditions everywhere except near the 


extreme tips. The aerodynamic inertia of the elliptic 
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wing is the same as if the sections were in two-dimen- 
sional flow except for the factor E. Hence: 


is ee = ra/E 

According to some approximate calculations (based 
on the lifting-line theory) made by the author,‘ the 
lift of the elliptic wing approaches its final value more 
rapidly than would be suggested by the two-dimensional 
theory. 

In the case of A = 3, or in the equivalent case of 
A = 6 with the wing vibrating symmetrically in torsion, 
the development of lift appears to be nearly instan- 
taneous. 

APPENDIX 


The problem of the fluid motion produced by trans- 
lation of a solid ellipsoid was first solved by Green in 
an investigation of the vibration of pendulums. For- 
mulas given for this problem in textbooks on hydro- 
dynamics become indeterminate when applied to the 
case of an elliptic disc. The following short discussion 
is therefore presented to show the application to the 
present problem. 

As explained in ref. 5, the surface potential of an 
ellipsoid can be given by a very simple formula. For 
motion along a principal axis, the potentia’ at any 
point on the surface is proportional to the coordinate 
of the point in the direction of motion. Thus, if the 
ellipsoid with semiaxes a > b > c along x, y, and gz, 
respectively, is moving with unit velocity in the direc- 
tion of z, the surface potential is simply 


¢ = C2 (3) 


The equipotential lines are the similar ellipses formed 
by the intersection of the ellipsoidal surface with a series 
of equidistant parallel planes perpendicular to z. The 
constant C depends on the axis ratio and its evaluation 
involves a special class of transcendental function 
The solution for the 


known as ‘“‘Green’s Integrals.” 


surface potential appears in the form: 
¢ = [y,/(2 — v,)]2 (4) 


where 
dx 


‘ = ave f — eames, —= = 
0 (c? + X) ya" + \)(b? + X) 


(c? + X) 


is Green’s integral. 
The reduction of these integrals to the standard 
Following equation 


elliptic functions is given in ref. 5. 
= (), in order to re- 


(6.1) of ref. 5 and substituting » 
strict the solution to the surface of the ellipsoid, will 


7 2 -2 = 

abi a* — ¢ bes r 

— — ———_ | y —c- | 
42 — ¢2)3 b? — c* Lac 


give 


aa 
Jo 


= yi 


where E is the complete elliptic integral with the 


modulus k = V (a? — 6*)/a®. The integral E is equal 
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to the perimeter of a quadrant of the ellipse ab divided 
by the semiaxis a.° 
Y—> 2 [1 — (c/b)E] 


Asc—> 0 


Since Eq. (4) becomes indeterminate ( —> 0 and 
Y, —> 2), it is necessary to express the solution in 


terms of x and y, which are related to z through the 
equation of the ellipsoidal surface: 


x? ¥* 2° , 
— + — + — as (5) 
a b? c* 


or 


aie elt 
Z=¢ }] = -—— — (6) 
\ a? 5? : 


Substitution for y, and z in Eq. (4) gives 


, ae a 


a ile -_ 
EW @ & ( 


or 
(x?/a*) + (y°/b*) + [¢°/0/E)"]}=1 (8) 


Hence the distribution of the surface potential over 
the disc may be represented by the ordinates of a cir- 
cumscribed ellipsoid having the vertical axis 2b/F. 
For infinite axis ratio, EK = 1 and the chordwise cross- 
sections of the potential distribution are circles of 
radius b. 

In order to illustrate the analogy to two-dimensional 
flow, it is convenient to introduce the angle @ defined, 
at a particular value of x, by 

cos 6 = y/y, (9) 


where y, = 6 /1 — (x?/a?) is the ordinate of the edge 


of the disc. Then, from Eq. (8), 


@ = (9,/E) sin 0 (10) 
which is the potential function of the two-dimensional 
case except for the factor 1/H. It follows then that 
the edge velocity is also reduced from that in two- 
dimensional flow by the factor 1/E. 
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On Propeller-Tip Interference Due to the 
Proximity of a Fuselage 


ALBERT GAIL* AND HO-SHEN LU! 
University of Minnesota 


SUMMARY ° 


The periodic changes of the air forces acting on a propeller blade 
element that passes by the flanks of a fuselage nose are theoreti 
cally investigated. This forcing function known to excite vibra- 
tions of propellers and airplane structure is harmonically ana- 
lyzed. The strength of the fundamental harmonic and the rela- 
tive strengths of the higher harmonics are found to be greatly 
dependent upon the location of the propeller plane downstream 
from the fuselage nose and upon the ratio of fuselage diameter to 
propeller diameter. 

The farther downstream the propeller is located behind the 
fuselage tip and the greater the fuselage diameter compared to 
the propeller diameter, the smaller is the intensity of the funda- 
mental harmonic and the smaller is the relative importance of 
the higher harmonics of these propeller-tip excitations. 

The clearance between propeller-tip and fuselage flank is found 
to have comparatively little beneficial effect within the usual 
limits at the designer’s disposal. 

The simplifying assumptions of this analysis restrict the 
validity of the results to fuselage shapes corresponding to a 
Rankine point-source half-body in flow of axial symmetry and 
propeller locations one fuselage diameter or more downstream 
from the fuselage bow. 

Experimental verification of this analysis and an investigation 
of oblique flow about the fuselage are desirable. 


INTRODUCTION 


O* MULTI-ENGINE airplanes having propeller blades 
passing the fuselage flanks (Fig. 1) strong vibra- 
tions have been observed, the cause of which could be 
traced to aerodynamic interferences between propeller- 
tips and fuselage. These interferences are the object 
of this analysis. 

The disturbances of the velocity field due to the pres- 
ence of the fuselage change the lift and drag of the 
propeller blade element whenever it passes through the 
vicinity of the fuselage. 

The absolute magnitude of the periodically changing 
forces acting on the blade will not be determined in this 
paper pending further development of the lift theory for 
finite wings and propellers in unsteady flow. How- 
ever, the relative magnitudes of these periodically 
changing air forces acting on the blade are of primary 
concern to the designer who wants to know how to de- 
crease these interference impulses rather than how 
strong they are. Moreover, he wonders how important 
the higher harmonics of these excitations may be since 
he realizes that the changes of the air forces on the blade 
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Fic. 1 propeller-fuselage combination 


investigated. 


are obviously not sinusoidally varying during each 
revolution of the propeller. 

This investigation is therefore concerned with the 
relative magnitude of the interference excitations and 
with the relative importance of their higher harmonics. 


SURROUNDING THE FUSELAGE 
NOSE 


THE VELOCITY FIELD 


The fuselage nose will be represented by a Rankine 
solid containing a single point source.{ The assump- 
tions implied in this approach are obviously: axial 
symmetry of flow and fuselage and a fuselage contour 
corresponding to that of the fluid surface separating the 


source flow from the superposed rectilinear flow (Fig. 


\ 


5). The velocity potential of this flow is 
¢g = Ux - (Q/4av/x? + r? (1) 

where x, 7 = cylindrical coordinates. 
X-axis = axis of symmetry with point source 


as origin. 
source intensity. 
velocity of the rectilinear superposed 

flow, 1.e., airspeed with which the 

airplane flies. 

Expressing the source intensity Q in terms of the 
diameter 2) which the fuselage asymptotically ap- 
proaches (Q = mpo?U») and substituting in Eq. (1) one 
obtains 


Un — (Uopo?/4V x? + r?) 


g= 


(la) 


t Lock, C. N. H. (R. & M. 1239) and Wieselsberger, C. 
(Z.F.M., September 28, 1932, page 533) have previously in- 
vestigated the interference between a propeller and is nacelle 
employing sources and sinks. 
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Hence the axial or x-component of the velocity is 








-¥.U, User = Ur+ du (2) 


Uu _— a 
ox 4(x? + r2)'/? 


where Aw is the increase of the x-component of the ve- 
locity due to the presence of the source flow (fuselage 
nose). 

The radial component of the velocity is 


Ov sie Uopo°r (3) 
or A(x? + r’) 


v, is so small as to be negligible for x > po. Moreover, 
the influence of v, upon the interference impulses is not 
investigated here. This influence is believed to be 
negligible for x > 0. The neglecting of v,, however, 
limits the scope of this analysis to locations of the pro- 
peller plane well aft of the source, 7.e., not farther for- 
ward than about one fuselage diameter downstream 


from the fuselage tip. 


THE PERIODIC CHANGES OF THRUST AND TORQUE 
ACTING ON A PROPELLER BLADE ELEMENT 


The change due to Aw of the flow at any blade ele- 
ment at a distance, a, from the propeller axis is appar- 
ent in the vector diagram Fig. 2. 


Au 


Uo 








- 


WA 


The velocity vectors on a blade 
element. 


Fic. 2 


The decrease of angle of attack due to Avw is 
Aa = tan! [Auwa/(Uo? + AuUy + w*a?)] (4) 


or for Au small as compared to U) as well as to wa, Aa 
is approximately 


Aa = waAu/(U)’ + w*a*) (4a) 


Furthermore the decreases in the lift and drag of the 
element are proportional to Aa, at least for the small 
values of Aa which occur in this problem. Thence the 
decreases in the thrust A7 and torque AQ acting on the 
blade element will be proportional to Aa as well as to 
Au 


AT = K,Aa = KyAu 


(5) 


AQ = K,Aa = Ko Au 


These con- 


where K;, Ko, K7, and Kg, are constants. 


THE AERONAUTICAL 





SCIENCES 


stants are influenced by propeller shape and operation 
conditions. They are particularly affected by the 
time rate of change of circulation. This analysis does 
not attempt to determine these constants and rather 
confines itself to determining the relative magnitudes 
of the interference excitations AT and AQ as affected 
by the relative dimensions of the propeller-fuselage 
combination. 

Au as encountered by the blade element of the 
propeller rotating with the angular velocity w is obvi- 
ously a function of wt, the angular displacement of the 
propeller about its axis. Fig. 3 shows a front and a 
plan view of propeller and fuselage illustrating the 
dimensions involved. 






= ae 


1 a / 
K 

i\ 8 L-acosat'! 
ie 
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Front and plan view of 


Fic. 3. 
propeller and fuselage nose illustrating 
the symbols employed. 


From Fig. 3 one reads: 


y* a® sin? (wt) + (L — a cos (wt))? 
(6) 
= L? + a*? — 2La cos (at) 


where L = distance between X-axis and propeller 


axis. 
L=«e&+p,+R. 
€, = propeller tip-fuselage clearance. 
p, = fuselage diameter in propeller plane. 
R = tip radius of the propeller. 


Substituting Eq. (6) in Eq. (2) one obtains Au as a 
periodic function of wf: 
Uopox i 
5) 
4[x? + L? + a? — 2La cos (wt)] ’ 
Thus the periodic changes of thrust and torque acting 
on the blade element 


AT KU ( Au/ Us) 


(8) 


AQ = KoUo(Au/Us) 
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These periodic forces may excite vibrations of the 
propeller, the engine, and/or the airplane structure. 
Since AT = fi(wt) and AQ = fo(wt) are non-harmonic 
(Eqs. (7) and (8) and Fig. 4) a number of modes of the 








os Tw 
THE V&LOCITY CHANGE Su/U, 

AS ENCOUNTEFED By THE 

BLAOE ELEMENT. 


The velocity change Au/ Uo as encountered by 
the blade element. 
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responding system may be excited simultaneously. 
These forcing functions of the propeller-tip interference 
will, therefore, be subjected to a harmonic analysis. 


NON-DIMENSIONAL PARAMETERS 


Before entering into the Fourier expansion one may 
express the geometric dimensions in units of some char- 


acteristic length, say R. Thus define 


po = yoR = half the greatest fuselage diameter. 

py = y,R = half the local fuselage diameter at a 
distance x downstream from the 
point source. 

x = &R = location downstream from the point 
source. 

e, = 8,R = propeller clearance = L — (p, + R). 

a = ¢R = distance of the blade element from the 


propeller axis. 
Eq. (6) may now take the dimensionless form 


Au ae vo 


+ 1 x Yx)? i ~ 
U le+atat wets 


2(1 + B, + vz)f cos (at) | (9) 


Let 
A=#+(1+8+ +2 
B = 1+ 8, + wt 
C = y7k/4 


Then Eq. (9) simplifies to 


Au/U,) = C(A — B cos wt) ~*/? (10) 


The fuselage radius p, at any section x may be de- 
rived as follows: 

Stokes stream function, y, for a point source in recti- 
linear flow is in cylindrical coordinates x and r 








y = l yy: — @ 7 od 
2 4r Vxt+ 7? 


Bo . see 
—=- l oe. = po’ l (2 + <a — ) 
2 4 V x? + 1? 


The fuselage corresponds to the streamline y = 0, 
whence, setting ry = p, to indicate the radius of the 
fuselage at any x, the shape of the fuselage is given by 


i _ x 
Ss = 2 ai(2 + = =) 
2 Ver4+ 72 


This equation may be solved for p, with the result 








=~ ya - CSE +] 


Po 


(11) 


The negative sign is physically meaningless as a nu- 
merical example would reveal. 

The contour of the fuselage nose (Rankine point- 
source half-body) as determined by Eq. (11) is plotted 
in Fig. 5. ; 
FORCING 


ANALYSIS OF THE PERIODIC 


FUNCTION 


HARMONIC 


The periodic changes of thrust and torque due to the 
proximity of the fuselage nose are found to be func- 
tions of wt and have the form 


f(wt) = G(A — B cos wt)~*” (12) 


where G, A, and B are constants for any particular 
airplane flying under a particular condition of operation. 
Eq. (12) may be transformed into 


f(wt) = H(1 + a? — 2a cos wt)” (13) 
having substituted 
a = (A/B)[1 — V1 — (B/A)?] (14) 
and 
H = G[(1 + a?)/A]”” 
The function (1 + a? — 2a cos wt)~* in which 
0 < a < 1 may be expanded* into the series 
(1 + a? — 2a cos wt)? = bo/2 + Bb; cos wt + 
be cos 2wt + b3 cos 38wt + ..... + 5, cos nwt = 
1 n= @ i 
52 b, cos mwt (15) 


* The authors are indebted to Dr. Isaac Opatowski for bring- 
ing to their attention the following references and to Dr. Dun- 
ham Jackson for mathematical advice. 

Tisserand, F., Traité de Mécanique Céleste, Vol. I, pages 273- 
277, Gauthier Villars, Paris, 1889. 

Burkhardt, H., Trigonometrische Reihen und Integrale, 
Enzyklopadie der mathematischen Wissenschaften., Vol. II, 1. 

von Zeipel, H., Entwicklung der Stérungsfunktion, Enzy- 
klopaidie der mathematischen Wissenschaften., Vol. IV, 2. 

The authors will be glad to supply the mathematical details 
of the expansion to any reader who may be interested. 
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Fic. 5. The contour of a fuselage nose corre- 


sponding to a Rankine point-source half-body. 


The first two coefficients by and 5; can be obtained by 
the following formulas. 


bo = (1 — at) (1 + a’) F(a "| ~ 
T 2 

b; = 

where 


4/2 ‘wt 
F(a. *) = f emer nig = a 
2 0 V1 — a? sin? (wf) 
complete elliptic integral of the first kind. 


- a/2 : a 
F(a *) -f V1 — a’ sin? (wt)d(wt) = 


complete elliptic integral of the second kind. 








For 0 < a < 1 the higher harmonics are 


be 2( 0 -+ “Yb — 3h 
a 


bs = C + Vi ™ : by 


3 a 


(18) 








6 = 
vo a oO 
or 
y. 9 
2n — 1 a 2n — 1 


Eq. (18) determines the coefficient of any higher 
harmonic term when the coefficient of the previous 
two successive harmonics is known. 

The coefficients }; have been calculated for0 < a < 1 
and are plotted in Fig. 6. 

The coefficients of the higher harmonics were deter- 
mined by means of Eq. (18) and are expressed as multi- 
ples of 5; in Fig. 7. 

Although Eq. (15) and thence }, de, 
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| 
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bn/b, = f(a). (Circled notations refer to b, 
and b; as used in text.) 


not express the complete forcing function it becomes 
clear from Figs. 6 and 7 that a is a major parameter 
determining the magnitude of the fundamental and 
the relative strength of the higher harmonics. 

a itself depends upon the geometric dimensions of the 
propeller—fuselage combination, 7.e., yo, x/po, §, Bx. 
(Eq. (14).) 

Fig. 8 shows a in a three-dimensional diagram as a 
function of yo, the diameter ratio of fuselage and pro- 
peller, and x/po, the aft location of the propeller plane 
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downstream from the point source. It is noteworthy 
that the aft location x/pp of the propeller plane and the 
diameter ratio of the fuselage and propeller have the 
greatest effect on a (Fig. 8) and on the magnitude of the 
harmonics (Figs. 6 and 7). The location of the pro- 
peller section, § = a/R, has a comparatively lesser 
influence. This holds true for the same reason for the 
propeller-tip clearance. 

It must be remembered that Eq. (15) is not the com- 
plete forcing function. Expanding AT and AQ ina 
Fourier series 


AT = K,Uo 0k (it2)" (* + bd; cos wt + 
4 A 2 
be cos 2wt + ..... ) (19) 
AQ = KgUo mai(2 + ot) at + b, cos wt + ..... ) 
4 A 2 
defining 
b, = me[ it fy" b 
n 4 A n 
AT = KyUo[(bo/2) + bi cos wt + be cos 2wt + 
bs cos 3uf + ..... (20) 
AQ = KoUo[(bo/2) + bi cos wt + ..... ] 


Thus the coefficient b; is a measure of the strength 
of the first harmonic of the tip interference excitations 
so long as the operation conditions, 7.e., Uo, r.p.m. and 
propeller pitch, are kept constant. Under the same 
conditions b,,/b; is a measure of the relative strength of 
the higher harmonics of these excitations. 

Fig. 9 shows b; vs. x/po and yo for the propeller-tip 
section (¢ = 1.0) and the 0.75R station (¢ = 0.75). 





ve 


fe 
© = f(¥>%) | A? 


Fic. 9. (The circled notation ),; refers to 


b; as used in text.) 


— 


b,,/b; is shown in Fig. 7 as a function of a@ which in 
turn depends on x/po, Yo, ¢ and 8, as shown in Fig. 8. 


RESULTS AND CONCLUSIONS 


The periodic changes of thrust, A7, and torque, AQ, 


acting on a blade element which passes by the flank of 
a fuselage nose are expressed as functions of the veloc- 
ity increase Au due to the proximity of the fuselage 
(Eq. (8)). The relative velocity increase Au/Up is a 
periodic function of the angle of blade rotation wt (Eqs. 
(9) and (10) and Fig. 4) and also depends on the rela- 
tive dimensions of the propeller-fuselage combination. 

The periodic functions, AZ’ and AQ, are then ex- 
panded into Fourier series (Eq. (20)) the coefficients of 
which, b,, are found to depend upon the relative di- 
mensions of the propeller-fuselage combination, 7.e., 
vo, X/po, £, (8,) as shown in Figs. 7, 8, and 9. 

These figures indicate the following practical results: 

(1) The relative dimensions of the propeller-fuse- 
lage combination affect the strength of the fundamental 
harmonic b; as well as the relative strengths of the 
higher harmonics b,,/b; in the same sense. 

(2) The greater the fuselage diameter as compared 
to the propeller diameter, 7.e., the greater po/R = ‘Yo 
the smaller are b; and b,,/b:. This result is valid only 
if one compares the excitations of propellers having 
the same diameters, blade angle, blade shape, ‘and 
r.p.m. in combination with fuselages of various diame- 
ters. The interference excitations can obviously not 
continue to increase as the fuselage diameter decreases 
to zero; they do, however, increase for yo between 0.5 
and 1.5 within which limits conventional airplanes 
are believed to stay. The diameters of fuselage and 
propeller are determined by the designer on the basis 
of more important considerations than the propeller- 
tip excitations. It would also be out of the question 
to change either one of these diameters in order to 
remedy intolerable vibrations due to these interfer- 
ences. ‘Yo must therefore be considered an unchange- 
able constant for a given airplane. 

(3) As might have been expected, the intensity of 
the interference excitations diminishes the farther away 
from the fuselage the blade element passes. Hence, the 
blade tip sections (¢ = 1.0) experience stronger excita- 
tions than for instance those at 0.75R (¢ = 0.75). 
Likewise beneficial decreases result from increased 
clearances between propeller tips and fuselage flank, 
B, = €/R. The effects of changes of the propeller 
clearance are disappointingly small considering that the 
clearance can usually be changed by a few inches at 
the most, which is a very small fraction of the propeller 
radius. Practically no remedy can be expected for 
undesirable vibrations due to these tip interferences 
by increases of the propeller-tip clearance within such 
limits. 

The variation of the interference forces along the 
blade (AT and AQ = f(f)) is of interest for the analysis 
of vibrations of the propeller blades. 

(4) The intensity of the interference excitations b, 
and b,,/b; is greatly dependent upon how far aft (down- 
stream) of the fuselage tip the propeller plane is located. 
The farther aft of the fuselage tip the propeller plane is 
located, 1.e., the greater x/po, the less intense are the 
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interference excitations. This statement holds only for 
values: of x/po > 1, 7.e., for propellers about one or 
more fuselage diameters aft of the fuselage tip. In the 
immediate vicinity of the fuselage tip, 7.e, —1 < 
x/po < -+1, the assumptions of the derivation of the 
forcing function (v, = 0) do not hold. 

It is of interest to the designer that fuselages protrud- 
ing far ahead between the propellers produce much 
smaller propeller-tip interference impulses than fuse- 
lages that extend only about one fuselage diameter 
ahead of the propeller plane. This may not be a suf- 
ficient reason to decide the extent by which the fuselage 
nose is to protrude. This point of view, however, de- 
serves consideration as an argument in favor of a far- 
protruding fuselage nose. 

The influence of the relative locations of fuselage tip 
and propeller plane, x/po, upon the intensity of the 
interferences is here determined for one particular 
fuselage contour only (Fig. 5). It must be expected 
that differently shaped fuselages would behave differ- 
ently in this respect, because the origin of the inter- 
ference is the velocity change Aw as encountered by the 
blade element and as produced by the presence of the 
fuselage. The rate of change of Aw in the direction 
normal to the fuselage surface obviously affects the 
sharpness of the increase and decrease of AT and AQ = 
f(wt); the rate of change of Au must be the greater 
the smaller the radius of curvature of the fuselage 
contour (p, = f(x)). 
centrifugal forces on sharply curved streamlines call for 


This must be so since the greater 
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equally greater centripetal pressure gradients which in 
turn are accompanied by rapid decreases of velocity 
normal to the streamlines in accordance with Bernoul- 
lis’ equation. 


In the light of this interpretation the decrease in in- 
tensity of the interference impulses for increasing values 
of x/po is largely due to the decrease of curvature (in- 
crease of radius of curvature) of the fuselage contour 
downstream from the fuselage tip. This leads to the 
speculation that a locally restricted deformation of the 
fuselage flanks where the propeller tips pass by may 
control the inflow disturbances Au in such a manner as 
to decrease the interference excitations materially. An 
experimental investigation of the effectiveness of such a 
simple remedy would be most desirable. 


Furthermore, this analysis must be amplified by an 
investigation of oblique flow about the fuselage which is 
ever present on account of the induced velocity of the 
bound vortex of the wing and the inclination of the 
fuselage itself in many flight conditions. Such an 
investigation is believed to yield to theoretical treat- 
ment along the lines of von Karman’s analysis of the 
flow pattern around airships in oblique flow. (Ref.: 
von Karman, N.A.C.A. Technical Memorandum No. 
574, Calculation of Pressure Distribution on Airship 
Hulls.) 

Finally, it may be mentioned that the results of this 
paper are also applicable to the analogous interferences 
of ships’ propellers. 


An Approximate Method to Predict the Tran- 
sition or “Flare” Flight Path in the Take- 


off or Landing of an Airplane 


OMER WELLING* 


Civil Aeronautics Administration 


SUMMARY 


The nature of the transition phase of the take-off and landing 
of an airplane is briefly discussed in which it is shown that, in 
the usual case, the two maneuvers are virtually identical physi- 
cally, since in each the airplane may be considered to be a free 
body being simultaneously accelerated by two mutually per- 
pendicular forces between similar end conditions. It is also 
shown that the time histories of the two accelerating forces are 
in each case so interrelated that one determines the other and 
that therefore the pilot may control the time history of one of 
the forces directly, in which case that for the other is auto- 
matically determined; or he may indirectly control the time 
history of the other force through his direct control of the one, 
in which case the necessary time history for the directly con- 
trolled force is automatically determined; but he cannot exercise 
independent control over both forces in any case. 

Typical time histories of these forces have been obtained for 
one airplane in take-off and for two airplanes in landing by 
means of step-by-step integration of the summation force equa- 
tions which may be written. From these solutions it appears 
that the possible time histories which will satisfy a given set of 
end conditions lie within fairly narrow limits. Based upon a 
simple geometrical approximation to the form of these typical 
time histories, a method is presented by means of which to 
predict the flight path in terms of the initial and final values of 
flight path and climbing or sinking speeds and of the thrust 
and/or drag characteristics of the airplane. 


INTRODUCTION 


6 be problem of establishing the dimensions of the 
minimum size airport into or out of which an air- 
plane may safely be operated has received considerable 
attention in the air transport industry during recent 
years. This has proved to be a complex and difficult 
problem involving the performance of the airplane, the 
behavior of the pilot, the geographical and meteoro- 
logical conditions at and surrounding the airport site, 
and other related factors. The basic item of per- 
formance involved in the problem is the flight path 
described by the airplane during the process of taking 
off and landing. 

The complete take-off of an airplane may be con- 
sidered to take place in three phases, viz.: the ground 
run, during which the airplane is accelerated horizon- 
tally from a standing start to a take-off or “unstick”’ 
speed at the point of leaving the ground; the transition 
from level flight, at the instant of ‘“‘unstick,” to steady 
climbing flight during which the airplane is accelerated 
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upward from zero rate of climb to the rate which exists 
at the climbing flight path speed chosen and is also 
accelerated along the flight path from the ‘‘unstick”’ 
speed to the climbing speed; and finally, steady climb- 
ing flight during which the airplane is assumed to move 
with uniform rectilinear motion along a flight path 
making an angle with the horizontal. 

Similarly, the complete landing may be considered 
to involve the following three phases: the approach 
glide which may be assumed to be a steady condition 
consisting of uniform rectilinear motion along a flight 
path making an angle with the horizontal; the transi- 
tion from steady gliding flight to substantially level 
flight at the point of contact, during which the airplane 
is accelerated upward from the rate of descent in the 
approach glide to substantially zero at contact and is 
also decelerated along the flight path from the flight 
path speed during the approach glide to substantially 
the stalling speed at contact; and finally, the ground 
run during which the airplane is accelerated hori- 
zontally from the contact speed to zero. 

All of the phases identified above, except the transi- 
tions, have in the past received sufficient attention that 
fairly generally accepted methods exist by means of 
which to calculate the dimensions involved in each. 
Such is, however, not the case for the transitions, the 
problem of predicting the flight paths for which is 
considered here. 


DISCUSSION 


The following discussion of the nature of these 
transitions is divided into a separate discussion of that 
which occurs during take-off and that during landing 
because, although the two may be considered to be 
dynamically similar- (each involving substantially 
simultaneous accelerations in two mutually perpen- 
dicular directions), one is chronologically the reverse 
process of the other and the pilot technique involved 
in each is believed to be fundamentally different. 
There follows a list of symbols which are used: 


A = aconstant 

B = aconstant 

C = rate of climb (ft. per sec.), positive upward; coefficient 
D = drag (Ibs.) 

E = energy (Ib., ft.) 

F = force (Ibs.), positive upward and flightwise 

G = velocity with respect to and parallel to the ground (ft 


per sec.), positive flightwise 
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L = lift (Ibs.); (subscript)—lift 


= an exponent 
= distance parallel to the ground (ft.) 


T = thrust (Ibs.) 

V = airspeed (ft. per sec.), positive flightwise 
W = gross weight 

a = acceleration (ft. per sec.? 


c (subscript)—climb, contact 


e (subscript)—effective 


g = acceleration of gravity = 32.2 ft. per sec.?; (subscript)— 
glide 

h = height (ft.), positive if measured upward; (subscript) 
horizontal; i.e., along the flight path 

k (subscript)—kinetic 

1 (subscript)—landing 

o (subscript)—initial 


p (subscript )—potential 


t = time (sec.) 

v (subscript)—vertical 

w (subscript )—wind 

8 = flight path angle (radians) = C/G 


TAKE-OFF 


The take-off transition or “‘flare’’ may be executed in 
several ways. The airplane may be gradually accelera- 
ted upward from ‘“‘unstick’” until steady climb is 
established at the chosen climbing airspeed or the pilot 
may allow the airplane to continue to accelerate along 
the flight path without gaining height until some speed 
in excess of the “‘unstick’”’ speed is attained before 
beginning the upward acceleration which produces 
climb. This initial speed in excess of ‘‘unstick’’ may 
be less than the climbing speed or greater, in which 
latter case the early stage of the transition is a ‘‘zoom.” 
There is reason to believe, however, that for a given 
set of end conditions the differences in the overall dis- 
tances required due to possible differences in the inter- 
vening execution are negligible for all practical pur- 
poses. It also appears from an examination of a great 
number of photographic records of take-off flight paths 
that in the typical case the initial flight path speed is 
less than the climbing speed and the climb is con- 
tinuously accelerated from zero at the beginning to 
that which results from a chosen climbing speed and 
throttle setting. The pilot technique required for a 
given throttle setting is so to manipulate the elevator 
control that the climbing airspeed is attained but not 
exceeded. 

The flight path for such a typical transition may be 
represented by the diagram in Fig. 1, which also illus- 
trates the forces acting upon the airplane. Since in 
the usual case 8, is of the order 10° or less it is assumed 
that sin 8 = 6 and cos 6 = 1.00 and the following 
summation equations may be written with respect to 
the flight path: 


At the beginning of the transition: 


G =G,=Vo.i Ve 


Cc =0 
> =L-—-W=0 
%3=T-D= 


ny 


h 
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At any point during the transition: 
G =G=V¥j; VJ, 
C=C 
> =~L-We=F, 


s,=T-D-pWw=F, 


At the end of the transition: 
fae (2 


G =G=-Y, 2 
C =C, 
> =L-W=0 


T—-D-—-BW=0 


> 
“hk 


A study of the above equations leads to the following 
conclusions: 

(a) The pilot’s control of flight path speed is largely 
indirect since the immediate response of the airplane to 
elevator control is a change in angle of attack which 
directly affects F, but, aside from the comparatively 
negligible possible change in D, affects F, only through 
the change in the BW term in the >, equation due to 
the change in C produced by F,,. 

(b) It follows that in order that the flight path 
speed be continuously accelerated from V, to V,, it is 
necessary that the climb be accelerated continuously 
from 0 to C, and at such rate as to produce the desired 
flight path accelerations. 

(c) It follows further that F, must increase from 
zero initially to some finite value and return to zero at 
the end of the maneuver, while in the same time and 
as a result thereof, Ff, must be reduced from 7 — D 
initially to zero at V,. 

(d) A formal solution of the equations for the pur- 
pose of predicting the flight path appears, if not im- 
possible, then at least likely to be so involved as to be 
of no practical value. 

A solution to these equations has been obtained by 
means of step-by-step integration for a typical air 
carrier airplane for a transition assumed to take place 
between 1.15 V, and 1.33 V,, which was subsequently 
checked by photographic records of flight tests in 
which good agreement was obtained. This solution 
required a running assumption of C,, the correctness of 
which assumption could be judged only after the solu- 
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(Vertical scale is of the order 5 times the 


horizontal.) 








TRANSITION 


tion had proceeded far enough to make it apparent that 
the end conditions would or would not be satisfied. 
Experience with this solution indicated that the time 
history of F, which would satisfy the end conditions 
lay within a comparatively narrow band of values at 
any time, although no difficulty in executing the 
maneuver was experienced by the pilot. It follows 
from conclusions (a), (b) and (c) above that F,, must 
be equally narrowly defined. Typical time histories 
of these two forces are shown in Fig. 2. The flight 
path approximates quite closely to a parabola. 
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Fic. 2. Typical time histories of F, and F, during the 


take-off ‘‘Flare.’ 


LANDING 


There appears to be somewhat less likelihood of 
appreciable variation in the manner of executing the 
landing “‘flare’’ than is possible in the case of take-off 
primarily because the freedom of choice on the part of 
the pilot is almost entirely limited to the selection of 
the conditions from which it is begun. He may select 
the steady gliding speed, the rate of descent (if power 
is used), and the height above ground at which to 
begin the transition. Once these three selections have 
been made, however, the pilot has little choice but to 
reduce the sinking speed to zero at the ground, inter- 
rupt the maneuver, or ‘“‘crack up.’’ For reasons analo- 
gous to those discussed under ‘‘Take-off,’’ the time 
history of F, which will accomplish this is fairly nar- 
rowly defined and there will accompany this decelera- 
tion of the sinking speed an uncontrolled but necessary 
deceleration of the flight path speed. 

The critical consideration in the selection of the 
initial conditions is to insure that this deceleration of 
the flight path speed does not reduce the speed to the 
stalling speed before reaching the ground. For this 
reason the maneuver is ordinarily deliberately so exe- 
cuted as to arrive at the ground at a flight path air- 
speed in excess of the stalling speed. The airplane 
may then be set onto the ground immediately or al- 
lowed to float just off it until the stalling speed is 
reached when it will settle on. 

A typical landing transition flight path is illustrated 
by Fig. 3, which also shows the forces acting upon the 
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Fic.3. Typical landing transition flight path. (Verti- 
cal scale is of the order 5 times the horizontal.) 


airplane. Making the same assumptions as in the 
case of take-off, these equations may be written with 


respect to the flight path: 
At the beginning of the transition: 


G =G,=V,7V. 


& 


C = C, 
> =~L—-W=0 
> =T-D-—BpW 


At any point in the transition: 


G =G 
Calf 

> =L-We=F, 
r,=T—D-pw=-F, 


At the end of the transition: 


G =G, 
C =0 

> =L—-W=0 
%=T-D=fF, 


These equations are identical with those written for 
the take-off transition and indicate a relation between 
F, and F, of the same nature as that in take-off. 
They also appear equally difficult of formal solution. 
A solution of the equations has been obtained by means 
of step-by-step integration for the air carrier airplane 
investigated under take-off and also for a typical light 
airplane. These transitions took place between 1.18 
and 1.06 V, and zero thrust was assumed. The time 
histories of F, and F, obtained by this means have 
been plotted non-dimensionally in Fig. 4 in such fashion 
as to indicate the relative magnitudes of these forces 
for each case. It should be noted that the form of each 
of these curves for each case is practically identical. 
The flight path for each case also approximates closely 


to a parabola. 


APPROXIMATE METHOD TO PREDICT FLIGHT PATH 


The following approximate solution rests upon the 
assumption that the actual time histories of F, and F, 
can be adequately represented by the nearest triangle 
of the same approximate form and enclosing the same 


area. These assumed equivalent histories are indi- 
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Typical time histories of F, and F, during the 
landing ‘‘Flare.”’ 

cated by the dotted lines in Figs. 2 and 4. 

tion consists of the following three steps: 


The solu- 


(a) The establishment of a relation between /, and 
S, (see Figs. 1 and 3) in terms of C, and the 
initial and final values of G. 

(b) Derivation of the equation for the parabola to 
which the flight path approximates. 

(c) The development of equations by means of 
which to calculate # and S from the thrust 
and drag characteristics of the airplane. 

The complete solution is presented for take-off but 

is indicated only for landing in the interests of brevity. 


TAKE-OFF (SEE Fic. 1) 


Assume these effective constant forces: 


F,, producing AC = C, — 0, int, 


F,, producing AG = G, — G,, in ¢; 


Then: 


AC = C, = (F,,g/W)t., or: F,, = C.W/gt, 


and: 
G, — G = (Fig W)t., or: 


F,, = WG, — G.)/st, 


eh 


AG = 


Therefore: 
Fa] Pay - C., (G, ei G,) (a) 
If the actual forces were in fact constant, the follow- 
ing relations would exist: 


h, = AC /2a,, = WC2/2gF.,, 


c 


since a = 


F/M = F,./W 


S. = Gote + (AG /2a,, = Got, + [W(G, — Gp)?/2gF,,] 


~&) 
Se 


Now: 


i= S, (<. > =) =" nC + Ge 


»). = 
2S,/(Go + G) 


& 
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Therefore: 
IG, W(G, — G,)? 
sx [1 22] - Maa @ 
G, + G, 2gF., 
_WG-G)) G+G _ WG - 6) 
7 2¢F., G, — G, 2gF., 
Then: 
i . WC? 28F., Coty 
h, Se a ee a a — = aoe 
2gF,, W (G,? 7 i) (G,? TH G,?) F,, 
But: 
F,. Fes = C,/(G, — G,) 
Therefore: 
; ( G, — G, & 
h, S. _— on) — —- al a ’ _ (b) 
G,2 — G,? ha G, + G, : 
Also: 
kh, = (a,,t.”, ay = F, gt,’ 2W 
and: 


S. = Got, + (a,,t.?/2) = Got, + (F,,t.?/2W) 


to (b) above will now be 
Referring 


The equivalent relation 
established for the assumed actual forces. 
to Fig. 2, it may be seen that these latter forces may be 
expressed in terms of the effective forces thus: 


F, = (4F,,t/t,)t (¢ = Otoét = #,/2) 
= 4F, — (4F,,/t,)t  (t = t,/2 tof = t,) 
and: 
F, = 2F., — (2F,,/t)t (t=Oto =4,) 
Then: 


C= Sf adt = g/w/f F,dt 


| 4F,. 
= 2/W | —*idt = 
t, 
2) 


a constant (¢ = Otot = t,/2) 


2eF, 
ai 
Wt, 


or, since C = 0 when ¢ = O, the constant = 0, and 


C = (2gF,,/Wt,)t? (t = 0 tot = #,/2) 
= g/W / ‘4F, — oe aay me Oey _ Oe os 
7 : ‘ W Wt, 

a constant (¢ = 1,/2 to i = £4) 

Now: 
C = (gF,,/2W)t, when t = ¢,/2 
Therefore the constant = 
(gF,,/2W)t. + (gF,,/2W)t, — 
(2gF,/W)t. = — (gF,,/Wht., 








TRANSITION 


and: 

C = (gF,,/W) [4t — 2(#/t.) —t] (¢ = t,/2 tot = t,) 

Also: 

G= fa,dt = g/W J Fydt = (2gF,,/W)f'1 — (t/t,)dt 
= (2gF,,/W)[t — (#/2t,)] + G (t = Oto t = t,) 


be t./2 
ri Cdt = (gF,,/W) | fie t*dt + 
be 
f 4t — 2(t?/t,) — tat 
te/2 


= GF./W){| ( 3t,)0 +.a constant fj a4 


and: 


hk, = 


| 2 — (2/3¢,)t — tt +a constant ff } 
The first constant = 0 since when ¢ = 0, h, = 0. 
When ¢t = t¢,/2, however, 


h, = (gF.,/W)(2t*/3t,)1=1./2 


= (gF,,/W)(t.2/12) 


Therefore the second constant = 


(t,.2/12) — (é,27/2) + (¢,2/12) + (t,?/2) = (¢,?/6) 
and: 
oF, f ) te/2 9 t2] ) 
h, = —: 3 +| 2" tf —t, + “| 
W 1F ii 3t, 6 jier2 
oF, /t,? 2 t? tt, 
= U2 — -tf — iF + — —- — + 
VW (‘; t 3 6 2 
t,2 4 t2 ¢,? 
12 2 6 


hk. = gF.t.?/2W 


te le 
f Gdt = J G, + (2gF,,/W) [t—(1/2t,)e]dt 
0 


= Git, + (2gF.,/W)[(é.2/2) — (t.2/6) + a constant] 


while: 


S. = 


Since S = 0 when ¢ = 0, the constant = 0 


and: 
5, = Gt, + (2gF,,/3W)t.* 


A comparison of these expressions for h, and S, 
with those obtained above when F, is assumed constant 
will show that the two expressions for h, are identical 
while those for S, differ in that the second term under 
the assumed actual conditions is */; + 1/. = ‘*/; of 
that when F.. is assumed constant. It follows that, 


for the assumed actual conditions: 


WC2/2gF,, 


h. = 
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and: 


5 : W(G, — G,)? 
S; = Got, + 4/3 X i. a “o = 
2gF., 
owe . (ys 
e+ a 
3gF., 
It also follows that the average value of AG is */; of 
that when F,, is assumed constant, 7.e.: 


S, = 3S, 
ele AG GCG G C+ 
G, + 4/, — G, 4+ — ae _ of: ¢ 
~ vo vo 
Substituting this value into the expression for S,: 
3G, 2W(G, — G,)? 
S, — = 7 ~ S, . 
G, + 2G, 3gF., 
Then: 


; ( 3G, ) 2W(G, — G,)? 
BE gern ae 
G, T 2G, 3cF, 


(& + 2G, — *) ‘ 
G, + 2G, Bi 


 2WG — G)* | Go + 2G 


3gF., 2(G, — G,) 
W(G, — G,)(G, + 2G,) 
3gF,, 
ini WC,? 3gF. . 
h,/S. = — oo aaa eorie 
22F,. 2W(G, — G,)(G, + 2G, ) Fi, 
3C2F,, 
2(G. — G.)(G. + 2G.) F., 
But: 
F.,/ Fe, = Ce/(Ge — Go) 
Therefore: 
3C.2(G, — G,) 3c. 
he|'Sy = wee me ae (1) 
2C.(G, — G,)(G, + 2G,) G, + 2G, 


2. h = AS* (see Fig. 1) 


Assume: 
kh = AS” 
Then: 
dh/ds = NAS‘~! 
Now: 


when S = S,, h = h, and dh/ds = B,= C,/G; t.e.: 


Q 


NAS,*— = C,, 


€ 
whence: 


A = C/G.NS*"' 
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Substituting back into the original equation: 


h, = CS" jens" ~* = C,S,/G,N 


or: 
N = C,S,/G,h,, or, since h,/S, = 1.5C./(G, + 2G,) 
N = (C./G,)X [(G + 2G,)/1.5G,] = (G, + 2G,)/1.5G, 
and: 

Go + 2Ge - 
h/h, = (S/S) *°™ (2) 


3. Dimensions (see Fig. 1) 

At the beginning of the transition the total energy 
with respect to the ground is the kinetic energy of the 
motion, 7.¢.: 


E, = MG,?/2 
At the end of the transition the total energy is made 
up of a kinetic and a potential portion, thus: 


E, = (MG,2/2) + Wh, 


The total energy change is, therefore: 


AE = E, — E, = w | a +a,| 
2g 

It is apparent that this change in energy has been 
produced by the effective value of the total accelerating 
force, (T — D), acting through the distance traveled 
by the airplane along its flight path. Consistently 
with the assumption that cos 6 = 1.00, this distance 
along the flight path is assumed equal to S,. It follows 


that: 


caus, 2 
(T — D),S, = W | + ne | 
2g 


or: 





Ww . —@ 
S. = | & 7 Ge + ne | 


(T — D), 2g 
But: 

h, = [1.5C,/(G, + 2G,)]S,, from Eq. (1) 
Therefore: 


(T — D).S, = (W(G2— G,*)/2g) + 
[1.5WC,/(G, + 2G,)]S, 








L.5WC, 
—— ) i Rent etl. I , =, 7 2 = 2\ /é 
«7 L de G, 4 ce S. W (G G, )/ 2g 
E — D),(G, + 2G,) - LW | ; W(G,2 — G,?) 
Sees S, = ———— 
G, + 2G, 2g 
W G2 — G,*)(G, + 2G, 
“< \Go+2G) 


2g * (7 ~ ih. &. + 94) ~1409 


The simplest approximation to (TI — D), appears 
to be: ;, 


(T — D), = [((T — D), + (T — D),|/2 
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and the step-by-step integration indicates this to be a 
valid approximation since the distance history of(77—D) 
is practically a straight line. This appears likely to 
be generally true in view of the following: 

(a) The relation of T to V is approximately linear. 

(b) <A curve of L/D vs. Cy, is concave upward be- 
tween any probable values of C, likely to occur at the 
end conditions. 

(c) The lift, which is equal to W at either end of the 
maneuver, is equal to W + F, between these end con- 
ditions. Since D = L/(L/D), the effect of this upon 
D is of opposite sense to that of the shape of the L/D 
vs. C; curve mentioned under (b). 

All of the expressions necessary to predict the take- 
off transition flight path, which have been developed 
here, are summarized below: 





h,/S, = 1.5C,/(G, + 2G;,) (1) 
h/h, = (S/S), where: N = (G, + 2G,)/1.5G, (2) 
W G,? — G,*)(G, + 2G, 
S.=—* ( ee (3) 
2g (T — D)A(G, + 2G,) — 1.5WC, 
where: 
(T — D), = ((T — D), + (T — D),)/2 
LanpINc—(SEE Fics. 3 AND 4) 
i h,/S; 


By a process analogous to that which appears under 
“Take-off” above, it can be shown that, with the origin 
at the beginning of the flare: 


h,/S,; = —1.5C,/(2G, + G) (4) 
2. (ty — h) = A(S; — S)® 
Similarly, it can be shown that: 
N = (2G, + G)/1.5G, 


and: 





2Go + Gi 
h/y, = 1—- (1 = 3) 15 : 
3. Dimensions 
Proceeding as under ‘“Take-off:”’ 
E, = (WG,?/2g) + Wh; 
and: 


E, = WG;?/2g 


a 2 
w | eS = .. | 
2g 


Then: 


AE = 


(T — D),S; = 





TRANSITION 


But: 





h, = [—1.5C,/(2G, + G)]S), from Eq. (4) 
Therefore: 
* _ WG?-G,?) Lae Cf 
t-Rhe +e Es 
| (T — D).(2G, + Gi) — sea 5 — WG? — G,*) 
2G, + G a 2¢ 


W (G;?? — G,*)(2G, + G)) 
2 l g)(2G, + Gi) _ (6) 


2g “* (T — D),(2G, + G) — 1.5WC, 





For reasons given under ‘Take-off,’ the average of 
initial and final values appears a good approximation 
for (T — D),. 

Summarizing, to predict the landing transition flight 
path: 

h,/ S; = —1.56,, (2G, + G;) (4) 
h/h, = 1 — (1 — S/S,)%, where: 
N = (2G, + G)/1.5G, (5) 
(G;? ~~, G,*) (2G, a G;) 


2. 
X(T —D).0G,+G)-15wc, 


— 





; V 
S= > 


FLIGHT 
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CONCLUDING REMARKS 


It may be noted that the method neglects such 
secondary effects as the reduction in induced drag due 
to proximity to the ground or the gradient with respect 
to height above the ground which is usually present in 
the wind velocity. Neither of these effects appears 
likely to make significant difference in the typical form 
of the time histories of the accelerating forces, upon 
the assumed simplification of which the validity of the 
entire Allowance for the change in 
magnitude of the forces due to these effects may be made 
in the selection of the proper values of (TJ — D), and 
The accuracy of the 


method rests. 


G to be inserted in the equations. 
results of any practical application of the method, 
however, must depend upon the precision with which 
the pilot can control the climbing airspeed in take-off 
or the approach gliding airspeed in landing. For this 
reason it appears that all such secondary effects may 
be neglected without introducing significant error. 
This feature is also common to the generally accepted 
methods for calculating the flight paths for the other 
two phases of the complete take-off or landing. 





Airplane Vibration Tests as Related to the 
Flutter Problem 
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THE FLUTTER PROBLEM (GENERAL) 


fren FLUTTER problem has been treated analytically 
by several authors but an attempt will be made in 
this paper to summarize the general problem as one of 
the airplane designer and to demonstrate how the air- 
plane vibration test may be employed as a powerful 
tool in its solution. 

The designer’s problem is to determine in the early 
stages of an airplane design the minimum flight velocity 
at which aerodynamically energized oscillation of a 
component airfoil surface or combination of surfaces is 
possible and, to incorporate design features which will 
assure that this critical flutter velocity will exceed the 
maximum velocity at which it is intended to operate the 
airplane. In short he must know (1) what factors gov- 
erned by design affect the flutter characteristics, (2) 
how to measure and evaluate these factors and (3) how 
to determine from the measured values of these factors, 
with a reasonable degree of accuracy, the value of the 
critical flutter speed inherent with his design. It can- 
not be said, however, that there are as yet available com- 
pletely satisfactory routine methods applicable to cal- 
culating the flutter speed of all aircraft components. 
The airplane designer can, however, by applying basic 
principles control certain of these measured factors by 
design so that there is reasonable assurance that the 
flutter characteristics of the airplane are at least quali- 
tatively influenced in the proper direction. 


Nature of Parameters 


Analytical treatment, reinforced by controlled labo- 
ratory tests as well as correlation with design and 
service history data has definitely established certain 
variables as those determining flutter characteristics. 
These may be broadly separated into: (1) natural 
frequency of the airplane structural components in per- 
tinent modes of motion; (2) elastic characteristics of the 
subject structures including location of flexural axes 
and nodes and loops of vibratory motion; (3) mass dis- 
tribution characteristics of the structural compo- 
nents including center of gravity location, moment of 
inertia, mass density and controt surface dynamic bal- 
ance features; and (4) damping characteristics includ- 
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ing material damping and friction damping. Of these, 
only the mass characteristics are readily calculable di- 
rectly from design data, it being most expeditious to 
measure the remainder by laboratory tests. 


Calculation of Flutter Velocity 


The more difficult phase of the problem, namely that 
of establishing methods and procedures by which the 
measured parameters may be employed in arriving at a 
value of the flutter velocity applicable to any given de- 
sign, has been and is being studied intensively by re- 
search and testing personnel. There is every indication 
that, in the very near future, a designer possessing rea- 
sonably accurate values of the flutter parameters ap- 
plicable to his design will have available methods and 
procedures by which he either can accurately calculate 
the critical flutter velocity or he can determine at least 
whether or not the airplane will be aerodynamically 
stable throughout the intended range of operating 


speeds. 
Measurement of Parameters 


It is hence apparent that the measurement of the 
flutter parameters is an important and essential prob- 
lem. It is the primary purpose of this paper, then, to 
outline the technique of the laboratory vibration test 
and to describe in what ways it contributes to solving 
the problem through the measurement of necessary 
data. 


NATURAL FREQUENCY DETERMINATION 
Technique 


The two principal methods of natural frequency de- 
termination are those employing the measurement of 
amplitude of response to a periodic vibratory force, and 
the measurement of frequency of free vibration resulting 
from the application of a step function force. The com- 
mon phenomenon of response of an elastic member to a 
periodic applied force is the chief tool in the measure- 
ment of natural frequencies. For all intents and pur- 
poses, with small system damping such as is the case 
with aircraft structures, the true natural frequency is 
very close to the frequency of maximum amplitude re- 
sponse. The observation and measurement of this 
latter is most convenient and satisfactorily accurate 
for practically all aircraft vibration testing. Free vibra- 
tion of structural components is accomplished by the 
sudden release of applied static loads, the imposition of 
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impact loads, or the discontinuation of periodic applied 
forces. 


Oscillators 


A periodic vibratory force may be applied to a struc- 
ture by means of an oscillator or shaker. Practically 
all of the shakers used in aeronautical work are of the 
mechanical variety employing mass unbalance as the 
source of the applied force. Simple rotary unbalance, 
double rotary unbalance with elements phased to give 
periodic translational force or couple, and Scotch yoke 
unbalance are the common types. A fundamental 
characteristic of the mass unbalance shakers is that of 
producing a force of a magnitude proportional to the 
square of rotary frequency. This necessitates a large 
unbalance for use at low frequencies and vice versa. The 
unbalance used generally ranges from .065 inch Ibs. in a 
small shaker suitable for control surface excitation to 
from 2.0 to 4.0 inch lbs. for shakers suitable for ex- 
citing engine and fuselage response. The outstanding 
advantage of the mechanical mass unbalance shaker is 
its self-contained nature; usually it is attached directly 
to the body being excited and driven by an electric 
motor through a flexible drive shaft. No attachment 
to any rigid reference body is required. 

The elastic or so-called “rubber band type’’ shaker 
employs an elastic member, one end of which is at- 
tached to the structure being excited and the other to an 
eccentric or Scotch yoke motion. The application of 
this type is not as simple as for the mass unbalance type 
as it requires a stationary reference point but it has the 
capability of producing large forces at low frequencies. 

The frequency characteristics of shaker driving 
mechanisms are very important. Due to the rapid 
change in amplitude of, and hence also in power ab- 
sorbed by, a vibrating body near resonance, the oscilla- 
tor should have suitable speed stability characteristics 
and ample reserve torque and power to give complete 
control at frequencies on both sides of the peak. 


Amplitude Measurement 


The measurement of amplitude may be accomplished 
in any of a multitude of ways ranging from the sim- 
plicity of passing a paper chart under the point of a pen- 
cil attached to the subject structure to the relative 
complexity of the electromagnetic pick-up with its at- 
tendant power supply, integrating amplifier and re- 
cording oscillograph. The method best suited varies 
with the application. 

The measurement of angular motion is a more diffi- 
cult problem. It can be accomplished by linear meas- 
urements at at least two points disposed in a plane 
normal to the axis of rotation of the surface. In gen- 
eral, it is necessary to measure the instantaneous rela- 
tive phase of the two parallel amplitudes in order to 
quantitatively measure the actual angular motion, as 
this latter is usually accompanied by some translational 
motion. The method which best yields this data is the 


oscillographic recording of the simultaneous output of 
electric pick-ups properly located. By proper treatment 
(addition and subtraction) the translation and rotation 
may be separated. If, however, there is a nodal axis 
within the structure (that is having no translational 
motion) the angular motion may be obtained directly 
from one linear measurement. Since it is often the case 
that the nodal axis does not fall upon the body, as for 
the common problem of the measurement of motion of 
a control surface the hinge line of which is moving 
translationally, and since the use of electrical equip- 
ment if available is always relatively complicated, an- 
other method has been developed which yields the de- 
sired data in a simpler and yet satisfactory manner. 
A plane mirror is attached to the rotating surface in a 
location which preferably has a minimum of transla- 
tional motion. A light beam is directed from some sta- 
tionary point onto the mirror and the image reflected 
upon a fixed screen or ground glass plate bearing a scale 
and located at a convenient distance from the mirror. 
The amplitude of the light spot travel is measured ,on 
the screen, which amplitude, combined with the dis- 
tance of the screen from the mirror gives an exact quan- 
titative measure of the angular motion. 

Frequency measurement is most conveniently made 
with an oscillographic system recording the output of 
an electric pick-up. A well calibrated constant-speed 
motor driven photographic oscillograph is very satis- 
factory. A calibrated cathode ray oscillograph also 
may be employed in combination with electric pick- 
ups for frequency determination. Stroboscopic de- 
vices also are applicable. 


Phasing 

The phase relationship between the motions of dif- 
ferent portions of airplane structures is a very necessary 
type of data for the proper selection and identification of 
any given vibration mode. A very simple yet ade- 
quate device provides a means of differentiating be- 
tween “in phase’ and “‘out of phase’ motions. Metal- 
lic probe points, between which are wired in series a 
portable battery and lamp, are held adjacent to the two 
portions of the airplane structure it is desired to phase 
after electrical bonding between the same has been ac- 
complished. By allowing the probes to make contact 
(t.e., ‘‘tap’’) once per cycle the relative phase is indi- 
cated by whether or not the lamp flashes. 

Stroboscopic devices also may be employed. Also, a 
cathode ray oscillograph, operating upon the output 
of pick-ups upon the structure at the two points in 
question will give the phase. The cathode ray will in- 
dicate any relative phase as will the stroboscope when 
used with a common reference motion. 


WING FLUTTER MODE NATURAL FREQUENCY 
DETERMINATION 


Wing flutter is covered chiefly by. four modes, 
namely: flexural-torsional, flexural-aileron, torsional- 
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aileron and flexural-torsional-aileron. As is implied 
by their titles the natural frequency of three different 
systems is involved: the wing flexural or bending fre- 
quency, the wing torsional frequency and the fre- 
quency of the aileron control system. In cases where 
other elastic systems are included in the wing structure, 
as for example a flap or float installation, other degrees 
of freedom are correspondingly added and should be 
considered in combination with those basic degrees of 
freedom which are pertinent. 


Wing Bending 


The wing bending frequency of a cantilever mono- 
plane or the wing tip bending frequency of a braced 
monoplane or biplane is most simply obtained by at- 
taching a shaker at or near the wing tip, located in a 
chordwise position determined to be near the flexural 
axis or axis at which an applied force causes only flexure 
of the structure with no attendant twist. Since the 
measurement of natural frequency is not critically 
affected by small variations in this respect the approxi- 
mate figure of 33 per cent chord is found generally 
satisfactory. The most important flexural frequency is 
that of the fundamental symmetrical mode in which the 
full wing, as a beam, vibrates in a half wave. The mass 
of the shaker must necessarily be small compared to 
the mass of the wing so as not to appreciably reduce 
the resonant frequency below that which would result 
with no additional mass. This wing mode usually is 
readily detected by observation but a resonance curve 
plot of the amplitude of motion at some point aids in 
the accuracy of detection. Many analytical treatments 
of the flutter motion in which this mode is considered 
treat the wing as a cantilever fixed at the airplane 
center line. It is found, however, that the actual mo- 
tion possesses symmetrical node locations somewhat 
outboard of the center line with finite vertical motion 
existing in the fuselage in way of the wing. The loca- 
tion of these nodes of course affects the resonant fre- 
quency and hence if there are supporting forces pres- 
ent forcing nodes at points other than those applicable 
to the desired mode, the measured resonant frequency 
will differ from that strictly applicable. Fortunately, 
it has been found that in the case of a conventional 
landplane this is not a large source of error in that node 
locations observed with the airplane, suspended upon 
a relatively soft sling, thus simulating somewhat free 
suspension, do not differ significantly from those ob- 
tained with the airplane resting upon its landing gear. 
It is considered, however, that the use of a proper sling 
is the better arrangement especially for float seaplanes 
which can otherwise be supported only very rigidly at 
the center line by means of handling trucks. 

The flexural natural frequency also can be excited 
by means of applying a shaker to the fuselage and ob- 
serving or measuring the frequency of resonant tip mo- 
tion in a similar manner. In this method, a larger 
shaker must be employed, and it is usually convenient 


to attach the shaker at the engine. In this way, the 
fuselage mode accompanying the desired wing mode 
can be studied. It should be pointed out that as soon 
as other elastic structural members are excited in se- 
ries (so to speak) with the surface in question such as is 
the case when the wing is excited from the engine 
through the fuselage, the resulting motions become more 
complex and the possibility of selecting a fallacious 
natural frequency due to the resonant response of a con- 
nected structure becomes more imminent. 

Natural wing bending frequency also can be obtained 
by the frequency measurement of free vibration ex- 
cited by imposing an impact load at the tip, releasing a 
static load or stopping a shaker attached to the struc- 
ture. Higher frequency wing bending modes, both un- 
symmetrical and symmetrical, can be studied in a simi- 
lar manner by shaker excitation. 


Wing Torsion 


The most important natural frequency concerned in 
the wing flutter study is that of wing torsion. The sig- 
nificant mode is the fundamental in which chordwise 
wing sections twist about some axis of oscillation form- 
ing an approximate one-quarter wave outboard of a 
node near the airplane center line with the maximum 
amplitude occurring at the tips. This mode may be 
excited by means of an oscillator and the resonant fre- 
quency observed or selected from a resonance curve of 
angular motion; or the frequency of a free vibration 
resulting from the release of an applied force or couple 
may be obtained with frequency measuring equipment. 

When using the resonance method, the applied shaker 
excitation should be such as to produce as pure a tor- 
sional motion as is possible. For this reason it has been 
necessary to employ lirfear force vibrators located off- 
set from the axis of rotation far enough to produce as 
high a ratio of couple to bending shear as possible. This 
may be accomplished by attaching the shaker to the 
trailing edge structure or to a boom rigidly clamped 
or attached to the wing beam or fittings and extending 
forward or aft from the surface. If the shaker is at- 
tached to the trailing edge care must be exercised so 
as to insure that the rigidity of the structural carry 
through to the main torque carrying member is suffi- 
cient to preclude local resonance and a fallacious value 
for the observed natural frequency. Ifa torsional boom 
is used, it must be designed also to be sufficiently rigid 
to eliminate bending resonance and should at the same 
time be light, with its mass disposed so as to add a 
minimum of moment of inertia to the chordwise section. 

The torsional frequency of a braced biplane wing cel- 
lule may best be excited by applying a shaker to the 
interplane bracing so as to act in a fore and aft direc- 
tion and displaced vertically from the effective torsional 
spanwise axis. In this way, the vertical bending mode 
response is avoided, there being negligible fore and 
aft response due to the inherent high rigidity in this 
direction. Similarly, in applying shaker excitation to a 
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monoplane wing, if the boom employed extends vertically 
above or below the surface, the same advantages will re- 
sult. Coupling with the wing bending may also be 
eliminated by supporting the wing along the desired 
nodal axis thus forcing the desired mode. 

The natural bending and torsional frequencies are 
affected by the disposition of concentrated masses such 
as items of equipment, armament or landing gear es- 
pecially those located well outboard and either having 
their masses disposed longitudinally or their center of 
gravity positions displaced a considerable distance from 
the elastic axis. For this reason the flutter character- 
istics of an airplane may vary considerably with differ- 
ent conditions of loading. Likewise noticeable is the 
change of natural frequencies resulting in a convertible 
airplane with a change from landplane to seaplane type 
gear. Here the tip floats, located well outboard and 
below the wing, considerably depress the natural fre- 
quencies. For such a case the bending frequency may 
be obtained in the same manner as with no tip float, 
but it is easier to excite the torsional mode by applying 
a shaker to the float or float bracing structure so that 
the force is directed fore and aft. To check for the 
effect of the flexibility of the float bracing in both the 
fore and aft and transverse directions as related to 
observations of resonant frequencies of wing torsional 
and bending response, respectively, the wing should be 
supported at the float bracing attachment points, and 
the natural frequencies of the float determined relative 
to the wing so fixed. 


Aileron System 


The aileron system, consisting of control surfaces, 
horns, bell cranks, push rods, cables and cockpit con- 
trols is a complicated combination of various masses, 
flexibilities and kinematic relationships. The mode of 
motion of interest in the flutter consideration is that of 
angular motion of the surface about its hinge line when 
restrained elastically by its control system. The im- 
portant mode is a symmetrical one which is independ- 
ent of the restraint supplied at the control stick, this 
being a node, and in which each aileron vibrates as a 
part of the elastic system formed by the controls on its 
side of the airplane. The natural frequency may be 
determined by employing shakers or by studying free 
vibration. A small mechanical shaker attached to the 
trailing edge of one aileron will satisfactorily excite the 
desired mode. It is found that no significant difference 
in observed symmetrical mode resonant frequencies 
results with the control stick free, hand held, or braced 
rigidly to the fuselage structure. It is convenient, 
however, to facilitate the measurement of vibratory 
response amplitudes by eliminating the superposition 
of extraneous angular excursions of the system by se- 
curing it rigidly in the cockpit. 

The problem of the increase in moment of inertia of 
the surface due to the addition of the shaker mass be- 


comes proportionally more important in the case of the 
aileron inasmuch as, with indiscriminate use of heavy 
shakers on small control surfaces, the resonant fre- 
quency observed will be considerably below the proper 
system natural frequency. To correct for the addi- 
tional mass contributed by the shaker, the resonant 
frequency of the system may be determined with sev- 
eral different amounts of additional mass, the results 
plotted logarithmically and an extrapolation made for 
the known amount of the shaker mass. This method 
may be used also in lieu of knowledge of the surface 
moment of inertia. 

With the controls restrained, a lower frequency un- 
symmetrical mode of vibration is observed in which 
the surfaces oscillate out of phase (7.e., one up and the 
other down) with some motion occurring at the stick. 
This frequency usually varies with the restraint, how- 
ever, and is not considered as significant in the flutter 
problem. When flutter calculations are made for the 
mode involving unsymmetrical aileron motion, the 
conservative assumption of zero natural frequency is 
made. Measurements of angular surface motion are 
obtained easily by the previously described mirror- 
light beam method. 

Ordinarily, the aileron system natural frequency is 
considerably lower than that of wing torsion so that 
coupling between the two modes is negligible. The 
bending frequency, however, is more likely to be close 
enough that it may be necessary to support the wing 
against deflection so as to eliminate coupling and the 
consequent shift of natural frequency from uncoupled 
values. It is noted that the reflected light beam 
method of angular measurement, when properly em- 
ployed, differentiates between angular and linear dis- 
placements and is not influenced by the presence of the 
latter. 


TAIL FLUTTER MoOpDE NATURAL FREQUENCY DETER- 
MINATION 
Tail Surface Modes 


The tail flutter modes involve both the coupling 
of the fixed and movable control surfaces in different 
degrees of freedom and the coupling of the tail surface 
with the various modes of the fuselage structure. Inter- 
action of symmetrical bending of the horizontal fixed 
surfaces with oscillation of the elevator system can 
occur in much the same way as wing bending-aileron 
flutter. Also, horizontal tail surface modes comparable 
with the wing flexure-torsion and torsion-aileron are 
within possibility but do not represent highly probable 
conditions. Transverse bending of the fin coupled 
with oscillation of the rudder system is a further con- 
sideration. A determination of the natural frequencies 
of the stabilizer bending and torsion, of the fin in bend- 
ing, and of the elevator and rudder systems supplies 
such necessary frequency data for the investigation of 
these modes. 
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Stabilizer 


The bending frequency of the stabilizer is most easily 
obtained by mounting a shaker on or near the tip and 
observing resonance. As with the wing, care must be 
taken in locating the shaker so that coupling of the de- 
sired mode with other modes is minimized. Also with 
small surfaces, in which the weight of the shaker is con- 
siderable relative to that of the surface itself, it may be 
necessary to apply it near the root, shaking the fuselage 
attachment points vertically so as to excite the proper 
flexural mode without the addition of excessive mass 
at the tips where it is very effective in reducing the 
resonant frequency. However, every resonant mode 
as shown by a resonance curve of amplitude measured 
at a point on one side of the fuselage must be com- 
pletely investigated by means of an amplitude and 
phasing survey to insure that the proper frequency is 
selected. By this means, other modes such as fuselage 
vertical bending, fuselage torsion, fin bending, unsym- 
metrical stabilizer bending, or motion resulting from 
local flexibility of fuselage attachment points can be 
identified. By supporting the conventional landplane 
upon its landing gear (including tail wheel) the vertical 
bending frequency of the fuselage is eliminated. If 
such local support of the tail is not feasible, the reso- 
nance frequencies obtained in the measurement of 
stabilizer motion must be carefully compared with 
those obtained in the measurement of fuselage motion 
at the tail in order that the proper frequency resulting 
in maximum ratio of stabilizer tip to stabilizer root 
The stabilizer natural fre- 
observations 


motion may be selected. 
quency also may be determined 
or measurements made while supporting the airplane 
from a sling and applying a shaker at the fuselage loca- 
However, as pre- 


from 


tion corresponding to the engine. 
viously mentioned, the possibility of selecting a falla- 
cious natural frequency due to the resonant response 
of some structure located between the shaker and the 
tail is much more probable. 

The torsional natural frequency of the stabilizer, 
usually very high in the conventional tail with its con- 
ventionally low aspect ratio and short span, is accord- 
ingly more difficult to excite and detect by shaker ex- 
citation. The methods of applying shakers offset from 
the nodal axis either by attaching directly to the surface 
or to suitable booms is applicable. 

Fin 

The fin bending frequency is similarly obtained by 
attaching a shaker to the fin tip so as to produce a 
force in a transverse direction and observing or meas- 
uring fin transverse amplitudes. In a manner similar 
to the technique described for the determination of the 
stabilizer bending frequencies, the same precautions 
must be taken of identifying the exact nature of all 
modes of resonant transverse response of the fin so as 
to differentiate the fin bending from the other related 
modes which also result in such motion. 
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Elevator and Rudder Systems 


The natural frequency of angular oscillation of the 
elevator or rudder about its hinge line as elastically 
restrained by the control system is the important con- 
sideration. This is obtained by attaching a shaker 
offset from the hinge line, preferably near the horn 
attachment points so as to eliminate torsional oscilla- 
tion of the surface itself. 

The elevator and rudder systems differ from the ai- 
leron system in that the elasticity and masses are not 
symmetrically disposed about the cockpit controls so as 
to cause the natural frequency of surface motion due to 
the elasticity of the system in the fundamental mode 
to be necessarily independent of the external restraint 
provided at the cockpit controls. Depending upon the 
kinematic and mass disposition characteristics of the 
control system itself, the actual method of restraint 
at the cockpit controls may or may not affect the value 
of the resonant frequency. The method of restraint 
most applicable is that which would be provided nor- 
mally by the pilot. In any event the frequency goes 
up with the rigidity of restraint and the upper limit 
is represented by the condition wherein the cockpit 
controls are rigidly braced to the fuselage structure at 
the point where the control force is normally applied 
by the pilot. Conversely, the case of minimum re- 
straint is that obtained with free controls. The possible 
range of natural frequencies will hence lie between the 
limiting values obtained by these test conditions. 


Tail Surface-Fuselage Modes 


Those possible tail flutter modes in which the motion 
of the fuselage constitutes one of the degrees of freedom 
include: 

(a) Fuselage torsion-rudder system: in which dy- 
namic mass unbalance of the rudder surface referred 
to its own axis and the effective nodal axis of fuselage 
torsional motion, causes coupling of the two modes. 

(b) Fuselage vertical bending-elevator system: in 
which the vertical motion of the tail portion of the fuse- 
lage is coupled with oscillation of the elevator system 
by lack of dynamic balance of the latter. 

(c) Fuselage transverse bending-rudder 
which is strictly analogous, to the above mode with 
motion occurring normal to rather than in the plane 
of symmetry of the airplane. 

(d) Fuselage torsion-elevator torsion: in which tor- 
sional deflection of the elevator, unsymmetrical with 
respect to the airplane center line is induced from fuse- 
lage torsional motion through dynamic unbalance of 
each half of the elevator surface relative to the hinge 
and fuselage torsional axes. If there is a local flexibility 
of the stabilizer relative to the fuselage about a fore and 
aft axis through the attachments, it is considered that 
this mode may be substituted for the fuselage torsional 
mode in that it produces the same type of elevator 


system: 


hinge line motion. 
It is to be noted that all of these tail flutter modes 


Lone 
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have been described as having two degrees of freedom. 
However, just as with wing flutter, those modes with 
more degrees of freedom are the more general possibili- 
ties. Hence for example fuselage transverse bending- 
fuselage torsion-rudder system flutter may be a con- 
sideration. 


Fuselage Modes 


The natural frequencies of vertical and transverse 
flexure of the fuselage structure are those at which 
resonant response at the tail in the pertinent direction 
results from excitation of the fuselage causing it to 
vibrate as a whole. They can be interpreted also as 
those frequencies of excitation at which the tail loca- 
tion occurs at a loop of the fuselage elastic curve. These 
data are best obtained by applying a large shaker to the 
airplane (most conveniently at the engine) and measur- 
ing the amplitude response at the tail. The resulting 
resonance curves, after complete identification of the 
airplane mode occurring at each peak, will indicate the 
proper fuselage natural frequencies pertinent to tail 
flutter. When vibrating vertically, the amplitude 
should be measured vertically on the airplane center 
line at or near the main stabilizer beam. When ex- 
cited by a transverse shaker force, the amplitudes 
should be measured transversely at the main fin beam 
at a vertical position which has been determined by 
preliminary survey to be the node of any torsional mo- 
tion which may be induced. The airplane should pref- 
erably be freely supported in a sling for such vibration 
tests. To eliminate extraneous control surface mo- 
tion, the cockpit controls are best restrained rigidly. 

The natural frequency of the empennage in torsion as 
elastically restrained by the fuselage can most directly 
be obtained by applying a shaker on the tail offset from 
the nodal axis so as to exert a force either transversely 
at the fin tip or vertically at the stabilizer tip. That 
resonant frequency of fuselage torsional amplitude is 
selected which a knowledge of related tail surface natu- 
ral frequencies and a complete phasing identification 
of fin and stabilizer tip motions shows to be the appli- 
cable mode. The tail should be freely supported for such 
an investigation but it has been found in many cases 
that the resting of a landplane tail wheel upon the floor 
did not alter the value of the observed torsional natu- 
ral frequency. It has been noted also that in the case 
of some airplanes having tubular fuselage construction 
in which a vertical section is a tall rectangle and which 
have different transverse rigidities in the upper and 
lower panels, that more than one natural torsional mode 
results. 


Elevator Torsion 


In this mode the opposite sides of the elevator os- 
cillate in an approximate quarter wave unsymmetrically 
about the airplane centerline as a node. The elevator 
in a conventional system is connected by a continuous 
torque carrying member and the connection to the 


control system is made at or near the airplane center 
line where the node occurs. Hence in order to excite 
this mode, the shaker should be located offset from the 
hinge line near one tip, and the horn should be prefer- 
ably restrained. However, if the subject natural fre- 
quency sufficiently exceeds that of the symmetrical 
mode of the system, the error due to coupling with 
control system elasticity will be small. If the excited 
torsional twist of the elevator tip is too small for accu- 
rate observation, the use of the reflected light spot 
method of angular measurement is convenient. The 
unsymmetrical mode should not be confused with the 
symmetrical torsional mode in which the tips twist in 
phase accompanied by an opposite twist of the central 
portion. This occurs at a higher frequency. 


DETERMINATION OF ELASTIC CHARACTERISTICS 


Certain of the elastic characteristics of aircraft 
structural components which are necessary for use in 
the investigation of flutter speeds are measurable in the 
laboratory. This consists of completely determining 
the deflection or elastic deformation characteristics of 
the structure in the various modes which are concerned 
in possible flutter motions. Included is the deter- 
mination of the nature of the elastic curves resulting 
both at zero frequency of applied force or couple and at 
a frequency corresponding to the natural frequency. 
When the designer uses certain natural frequencies in 
his flutter calculations he should consider fully the 
manner in which they were obtained, whether they are 
for coupled or uncoupled modes, where the nodal axes 
are located, and how such locations compare with that 
of the flexural axis. 

The elastic or flexural axis has been defined as that 
axis about which the surface rotates under the applica- 
tion of a pure static couple. It is also defined as the 
chordwise position at which an applied static force 
will produce only pure translational displacement of the 
surface. The location of this axis may be determined 
by two methods—those implied by the alternate defini- 
tions. A pure couple may be applied near the tip sec- 
tion and measurement made of the locus of chordwise 
points which undergo no displacement. Alternately, 
a force may be applied in successive chordwise loca- 
tions at each wing section and that location observed 
at which the vertical displacement of leading and trail- 
ing edges are identical. The wing must be securely sup- 
ported and the measurements made accurately in 
order to produce useful data. Coincident with the 
former method, a measure of torsional stiffness useful 
in divergence speed calculations may be obtained from 
the known values of applied couple and measured values 
of torsional deflections. 

A very important use for the knowledge of nodal 
axis locations concerns all modes involving control 
surface flutter and that is its use as the basis for dy- 
namic balance requirements for the surfaces. An axis of 
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torsional fuselage motion observed in the laboratory is 
a much safer basis for dynamic balance calculations for 
the rudder than an arbitrary center line of thrust loca- 
tion or center line of supporting structure. Node loca- 
tions of fuselage bending observed in the laboratory 
will be the most logical and conservative bases for bal- 
ance calculations inasmuch as their selection will have 
taken into account the actual elastic characteristics of 
the given structure including the effect of cut outs. 

The actual flexural mode of a surface such as a wing 
or horizontal tail directly affects the natural frequency. 
If ‘‘unnatural”’ node locations are forced due to rigidity 
of supports, ‘‘unnatural’”’ natural frequencies result. 
If the node location is artificially forced at the wing root 
such as would be obtained by rigidly supporting the 
fuselage or by use of a jig, the natural frequency is 
lowered. Also, if nodes are artificially forced outboard 
of their natural location by rigid supports or by sup- 
port by landing gear running to outboard points on 
the wing, the natural frequency is likewise lowered. 


DETERMINATION OF MASS CHARACTERISTICS 


Of the pertinent mass characteristics involved in 
the flutter consideration, none are primarily determined 
by laboratory tests as analytical calculation ordinarily 
produces sufficiently accurate data. The function of 
the laboratory test in checking such calculated data, 
however, merits some mention. 

Center of gravity or inertia axis positions of fixed 
surfaces can be checked approximately by the location 
of the nodal axis in the torsional frequency determina- 
tions as previously mentioned. Movable control sur- 
faces can be removed and readily checked for c.g. posi- 
tion. By swinging on a calibrated torsional pendulum 
the moment of inertia may be determined. A method is 
also available for laboratory quantitative determination 
of product of inertia by forced high frequency oscillation 
of the surface about one axis and observing the location 
of nodes resulting in the motion induced about the other 


axis. 


DETERMINATION OF DAMPING CHARACTERISTICS 


Internal structural damping is an important factor 
influencing the flutter phenomenon. It is known that 
an increase in the damping characteristics of any elastic 
system concerned in flutter will raise the critical velocity 
at which flutter is encountered. In the absence of data 
concerning applicable values of damping characteristics 
it is conservative to calculate the flutter speed assuming 
zero damping. However, the difference in calculated 
speeds between an assumption of zero damping and one 
of a small amount of damping may in some cases be 


great. Hence, if zero damping leads to a flutter speed 


value not sufficiently high, it is reasonable that the 
designer is concerned with determining values of the 
damping characteristics which are applicable to his de- 
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sign and which will result in a reasonable yet conserva- 
tive value of flutter speed. 

It is important to note that in the analytical treat- 
ments of flutter the aerodynamic forces and moments 
upon the airfoil in its non-steady motion are com- 
pletely expressed so as to include the oscillatory com- 
ponent of motion, and hence the damping coefficient 
correctly applicable to the calculations must represent 
only the internal dissipative forces and should not in- 
clude any of the aerodynamic damping necessarily en- 
countered in the oscillation of an airfoil in the at- 
mosphere either at rest or at any forward velocity. 
Most of the laboratory methods of measuring damping 
unfortunately involve dynamic vibratory tests and 
hence include aerodynamic effects. 

Internal structural damping may be evidenced as (1) 
creep or plastic flow, (2) elastic hysteresis or deviation 
of a homogeneous material from linear elasticity, ex- 
emplified by the area enclosed by a stress-strain curve 
and (3) friction or relative motion of parts of a non- 
homogeneous structure usually occurring in bearings or 
mechanical joints. The inherent nature of aircraft 
structural materials eliminates the first from considera- 
tion. The effective combination of the latter two types 
may be evaluated by obtaining a load-deflection curve 
for the subject structure as it is taken slowly through a 
complete cycle of reversal of load. Because of the very 
low damping of a typical aircraft structure, extreme ac- 
curacy of measurement of load and deflection is essen- 
tial. The loading shouJd be such as to produce a de- 
flection curve approximating that of the mode under 
consideration. It is considered that concentrated loads 
produce a satisfactory approximation for the funda- 
mental modes. The damping coefficient may be evalu- 
ated from the energy represented by the area enclosed 
in the elastic hysteresis curve, the spring constant of 
deflection at the measured point and the maximum am- 
plitude of cyclic deflection. 

Damping may be determined from displacement- 
time curves of free vibration obtained by the sudden 
release from an initial deflection or the sudden removal 
of a vibratory shaker force occurring at a frequency 
near that of free vibration. Damping may also be de- 
termined in several ways from resonance curves of 
motion of the structure in the pertinent mode. The 
damping coefficient can be determined from certain 
well known relationships involving the phase angle be- 
tween vibratory force and actual displacement as meas- 
ured at frequencies bearing certain definite relation- 
ships to the resonant frequency. The ratio of the 
maximum response at resonance to the static deflection 
that would be obtained under an applied load equal to 
that exerted by the shaker at resonance is also an in- 
verse measure of the damping. 

It is acknowledged that some of these relationships 
indicated above apply only to systems having low orders 
of damping, but it is considered that aircraft structures 
are satisfactorily included within this classification. 
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INSTITUTE MEETINGS 


Wednesday, December 17, 1941 


Wright Brothers Lecture. ‘‘New Pathways in Aeronautical 
Theory,” by Prof. Richard V. Southwell, Professor of Engineering 
Science, Oxford University, England. At the Pupin Physics 
Laboratories, Columbia University, New York City; 8:30 p.m. 

Dinner to Wright Brothers Lecturer. At Columbia University 
Faculty Club, New York City; 6:30 p.m. 


Tuesday, January 27, 1942 


Honors Night Dinner. At the Hotel Waldorf-Astoria, New 


York City; 7:30 p.m. 


Wednesday-Friday, January 28-30, 1942 


Tenth Annual Meeting. Annual Business Meeting of Mem- 
bers, Wednesday, January 28, 5:00 p.m. Three days of technical 
sessions—morning, afternoon and evening—on Aerodynamics, 
Airplane Design, Air Transport, Materials, Meteorology, Physio- 
logical Problems, Power Plants, Radio and Instruments, Rotat- 
ing-Wing Aircraft, Structures. 

The dates and hours for the various technical sessions will be 
announced in the December issue of the Journal together with 
the titles of the papers to be presented. A large number of 
papers have already been received and submitted to the Chair- 
men of the sessions concerned, who will shortly make the final 
choice of the papers to be presented at their sessions. However, 
any members who have prepared technical papers or contemplate 
doing so are still urged to send their papers or outlines to the 
Institute. Even if it should be found impossible to include their 
papers in the appropriate session at the Annual Meeting, the 
session Chairmen and the Editorial Board would like to have 
these papers to be considered for presentation at meetings of 
Institute Sections or for publication in the Journal. 


New CorporaTE MEMBERS 


The Institute is pleased to announce the following new Corpo- 
rate Members: 

American Bosch Corporation, Springfield, Mass., manufacturers 
of magnetos for four-, seven-, nine-, fourteen- and eighteen-cyl- 
inder aircraft engines. The company also manufactures fuel 
injection equipment. 

Dowty Equipment Corporation, Long Island City, N. Y., manu- 
facturers of the ‘‘Live-line’’ automatic hydraulic pump, the 
levered suspension type landing gear strut manufactured under 
license of Dowty Equipment Ltd. of England, actuating cylinders, 
valves and allied equipment. The company is also engaged in 
the development of a new type landing gear for the Air Corps. 


» News or InstiruTE MEMBERS 


Glenn D. Angle, A.F.I.Ae.S., Editor of the aviation yearbook 
“‘Aerosphere”’ and consulting engineer on aircraft engines, has 
organized at Tecumseh, Michigan, the Angle Engineering Com- 
pany of which he is Vice-President and Chief Engineer. The 
company will be engaged in the design, development and manu- 
facture of internal combustion engines, principally for aircraft. 

Gerald H. Backer, M.I.Ae.S., has joined the faculty of the 
College of Engineering, University of Kentucky, as Assistant 
Professor of Aeronautical Engineering. He had previously been 


Instructor in Engineering at the Casey Jones School of Aero- 
nautics and Academy of Aeronautics since 1934. 

Edward E. Brush, M.I.Ae.S., has become Associate Professor 
and Head of the Department of Aeronautical Engineering at the 
University of Kansas. He had been a Senior Project Engineer 
for Lockheed Aircraft Corporation and before that an Instructor 
in Aeronautical Engineering at the University of Minnesota. 

Pilot Officer C. G. Bull, Technical Member, formerly Mechani- 
cal Engineer with the Royal Canadian Air Force, is now assigned 
to flight instruction on advanced aircraft in the R.C.A.F. 

Adolph Burstein, M.I.Ae.S., Project Engineer with Vultee 
Aircraft, Inc., has transferred from the Stinson Division in 
Michigan to the company’s main plant in California. 

Col. V. E. Clark, F.I.Ae.S., has been appointed a member of 
the National Panel of Arbitrators of the American Arbitration 
Association, according to announcement of its President, C. V. 
Whitney, M.I.Ae.S. This organization expedites the settlement 
of commercial and industrial controversies through arbitration. 
Colonel Clark, now technical consultant to Hughes Aircraft 
Company, has been prominently active in the aeronautical engi- 
neering field since 1915 when he was Chief Aeronautical Engineer 
for the U.S. Army Air Services. He is well known for his de- . 
velopment of the “Clark Y”’ series of airfoils as well as for more 
recent work on the design of many types of aircraft, including 
plastic molded airplanes. 

Herbert F. DeWeese, Technical Member, formerly Enginee1 
and Aircraft Instrument Instructor at the Lincoln Airplane & 
Flying School, is now Instrument Installation Design Engineer 
with the Wichita Division of Boeing Aircraft Co. 

Herbert M. Heuver, M.I.Ae.S., has joined the Experimental 
Engineering Section, Materiel Division of the U.S. Army Air 
Corps at Wright Field as an Aeronautical Engineer. Except 
for a three-year period in 1937-39, Mr. Heuver has been an 
Engineer with the University of Washington since 1933. He 
designed the wind tunnel and other aeronautical equipment 
there in collaboration with Prof. F. K. Kirsten, A.F.I.Ae.S., 
and was in charge of the research work conducted by the Aero- 
nautical Department. 

Evert C. Hokanson, M.I.Ae.S., has been appointed to repre- 
sent the Aviation Department of the Whiting Corporation, 
Harvey, Illinois, in their recently opened West Coast office. 
Mr. Hokanson was previously engaged in the development of 
aircraft performance standards and standard flight test pro- 
cedures for the Engineering Flight Test & Inspection Section of 
the C.A.A. He has had experience as an Aeronautical Engineer 
with United Air Lines, Hall Aluminum Aircraft Company and 
Boeing Aircraft Company. 

Roger F. Horton, Technical Member, formerly an engineer 
with the Acrotorque Company, has been Chief Inspector for 
Dowty Equipment Corp. since June of this year. 

D. C. Jackson, Jr., M.I.Ae.S., has taken a year’s leave of 
absence from the University of Notre Dame where he has been 
Dean of the College of Engineering and Professor of Mechanical 
Engineering since 1939. He will go on active duty as a Major 
in the U.S. Army and will be stationed at Frankford Arsenal in 
Philadelphia. Major Jackson was a Lieutenant in the U.S. 
Army during World WarI. He has devoted his entire time, since 
receiving his S.M. degree from M.I.T. in 1922, to engineering 
education. This included posts as Director of Lewis Institute 
in Chicago, Head of the Electrical Engineering Department at 
the University of Kansas, and Head of the Department of 
Mechanical and Electrical Engineering at the University of 
Louisville. 

Sgt. Pilot Desmond F. Macdonald, Technical Member, has 
completed his Student Pilot course in the Royal Canadian Air 
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Force and is now a Staff Pilot of the Second Bombing and 


Gunnery School. 

Capt. S. G. Nordlinger, Technical Member, recently promoted 
from rank of Ist Lieutenant, is with the U.S. Army Air Corps 
Materiel Division at Wright Field doing experimental work on 
aircraft engines and accessories. He was formerly Chief Test 
Engineer, Detroit Division, Continental Motors Corp. 

James D. Redding, A.F.I.Ae.S., Formerly Assistant Chief of 
the C.A.A. Aircraft Airworthiness Section, is Staff Representa- 
tive for Aeronautics with the Society of Automotive Engineers. 

George Tharratt, M.I.Ae.S., was recently appointed Chief 
Engineer of Adel Precision Products Corp. in Burbank. He 
brings to this new position twenty years of experience in aircraft 
design and production in England, Canada and the United 
States. For the past two years he was Supervisor of the Depart- 
ment of Production Engineering which he formed at the Douglas 
Aircraft Company. During his association with Douglas, which 
he joined in 1936 as an airplane designer, he also organized a 
section to design and create interior furnishings for commercial 
transport planes. Previous to this he was an engineer with 
Great Lakes Aircraft Co. and The Glenn L. Martin Co. 

Frederick R. Weymouth, A.F.I.Ae.S., has joined the Depart- 
ment of Design Engineering of Bell Aircraft Corporation. He 
was recently Assistant Professor of Aeronautical Engineering at 
New York University and for the ten years preceding had been 
in the Engineering Departments of Curtiss-Wright Corporation 
at their Buffalo and St. Louis plants. 

Eli A. Zeitlin, Technical Member, has left the staff of the 
N.A.C.A. at Langley Field, where he was a Junior Engineer, to 
Assistant Aeronautical Engineer in the 


take a position as 


Bureau of Aeronautics, Navy Department. 
SECTIONS 


Los Angeles. Wing Commander Michael Crosley, represent- 
ing the British Air Commission at the North American Aviation 
plant, was the guest speaker at a dinner meeting of the members 
in the Los Angeles region held on October 14. Wing Com- 
mander Crosley is now a Test Pilot for the Royal Air Force, but 
had previously been a Fighter Pilot and has 35 victories to his 
credit. He is an authority on the subject of war in the air and 
spoke to an interested audience of over 250 Institute members 


on his experiences. 


New York. On October 8, a dinner meeting was held at the 
Advertising Club in New York City. This first meeting of the 
season for members in the Greater New York area was attended 
by 160 members and guests. After dinner, Victor J. Skoglund, 
Project Engineer of Pratt & Whitney Aircraft, read his paper on 
“Teing of Carburetor Air Induction Systems of Airplanes and 
Engines,” This paper is pub- 
lished in full in the October issue of the Journal. T. P. Wright, 
Past President of the Institute and formerly Director of Engi- 
neering for the Curtiss-Wright Corporation, gave a talk on ‘‘The 
Aircraft Program of the OPM and Other National Defense 
Mr. Wright is now Assistant ‘Director of the Aircraft 
He gave an interesting account of the mutual 


which was illustrated with slides. 


Agencies.”’ 
Section of OPM. 
problems faced by industry and government in the National 
Defense effort and the plans and procedures for their solution. 
His talk was followed by a lively discussion period. 

The Committee in charge of the meeting was composed of 
Dr. R. P. Harrington, Polytechnic Institute of Brooklyn, Prof 
F. K. Teichmann, New York University, and Walter M. Hartung, 
Casey Jones School of Aeronautics. Dr. Harrington was Chair- 


man of the meeting. 
STUDENT BRANCHES 


New officers for this 
Samuel Burton, 


Aero Industries Technical Institute. 


Branch are: Edwin Gustafson, Chairman; 
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Vice-Chairman; Philip Hier, Secretary-Treasurer. At a recent 
meeting, William Douglas of the Douglas Aircraft Company 


spoke on plant experience. 


Agricultural and Mechanical College of Texas. The first 
meeting of the school year 1941-42 was held on September 23 at 
which Prof. Howard W. Barlow, Honorary Chairman, presided 
Officers were elected as follows: Roy Grobe, Chairman; R 
Skidmore, Vice-Chairman; S. Kaffer, Secretary; and John 
Pratt, Treasurer. On October 9, Edgar Lesher read the Student 
Branch Lecture on ‘‘Typical Research Problems of United Air 
Lines” by R. D. Kelly. 


University of Alabama. Officers for the 1941-42 school year 
are: Angelo Conti, Chairman; Andrew Kuzmik, Vice-Chairman; 
and George Alderdice, Secretary-Treasurer. 

Boeing School of Aeronautics. At a meeting held on Septem- 
ber 3, officers were elected as follows: Raymond Tuthill, Chair- 
Robert Crowther, Vice-Chairman; Stanley Soult, Secre- 
Dr. Maurice Garbell is Honorary Chairman. 


man; 
tary-Treasurer. 


Indiana Technical College. On October 1, at its first meeting, 
officers for this Student Branch were elected as follows: Alan 
Starr, Chairman; William Hazard, Vice-Chairman; Norman 
Cosner, Secretary-Treasurer. Prof. James Ross is Honorary 
Chairman. 


North Carolina State College. A smoker was held on Septem- 
ber 30, the first meeting of the school year, at which Prof. L. R. 
Parkinson, head of the Aeronautics Department, was guest 
speaker. Officers were elected as follows: W. S. Johnson, Jr., 
Chairman; J. R. Weaver, Jr., Vice-Chairman; A. Michaels, 
Secretary-Treasurer. 


University of Notre Dame. Officers for this Branch are: John 
L. Hoelscher, Chairman; William Ungeschick, Vice-Chairman; 
Robert S. Mullaney, Secretary-Treasurer. 


Rensselaer Polytechnic Institute. Capt. W. M. Masland, 
Atlantic Division, Pan American Airways, spoke informally at a 
meeting held on October 9. The meeting was attended by ap- 


proximately 200 members and guests. 
PERSONNEL OPPORTUNITIES 


The Personnel Bureau serves individual members, as well as 
organizations seeking to employ aeronautical specialists. Any 
member or organization may have requirements listed without 
charge. 

Wanted 


Layout Draftsmen. Nationally recognized aircraft plant east 
of the Mississippi—a volume manufacturer of planes even in 
desires experienced layout draftsmen. Please list 
Address reply to Box 142, 


normal times 
your experience and salary required. 
Institute of the Aeronautical Sciences. 
Two Project Engineers and Several Aircraft Layout Draftsmen. 
Middlewest. Permanent posts to suitable men. In reply please 
give academic training, previous experience and salary desired. 
Address reply to Box 148, Institute of the Aeronautical Sciences 


CHANGES IN MEMBERSHIP 


The following changes in the membership of the Institute have 
occurred since the publication of the previous issue of the Journal. 


ELECTED TO MEMBER GRADE 


Cone, John Carrol, Asst. Vice-Pres., Pan American Airways 


System, Inc. 
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Doane, Foster Baird, B.S. in Chem. Eng.; Vice-Pres. and 
Treas., Maganaflux Corp. 

Nelson, Donald Earl, M.S. in Ae.E.; Stress Analysis Engineer, 
Northrop Aircraft, Inc. 

Vanderlipp, Nathan Floyd, Chief Engineer, Fairchild Aircraft 
Ltd., Canada. 


TRANSFERRED TO MEMBER GrabDrE 


Smith, Appollo Milton Olin, M.S. in Ae.E.; Asst. to Chief 
Aerodynamicist, El Segundo Div., Douglas Aircraft Co. 


ELECTED TO INDUSTRIAL MEMBER GRADE 


Fowler, Robert Erwin, Public Relations Director, Air Asso- 
ciates, Inc. 

Litchfield, Burnham, Vice-Pres. and Gen. Mgr., Allied Aviation 
Corp. 

Potter, Charles Allison, Aviation Representative, The Texas 
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